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CHAPTER I 

EXECUTIVE SUMMARY 

This work was aimed at exploring means for utilizing computer technology in 
quality inspection and evaluation. Inspection of metallic welds was selected as the 
main application for this development and primary emphasis was placed on visual 
inspection, as opposed to other inspection methods, such as radiographic 
techniques. Emphasis was placed on methodologies with the potential for use in 
real-time quality control systems. 

Because quality evaluation is somewhat subjective, despite various efforts to 
classify discontinuities and standardize inspection methods, the task of using a 
computer for both inspection and evaluation was not trivial. The work started out 
with a review of the various inspection techniques that are used for quality control 
in welding. Among other observations from this review was the finding that most 
weld defects result in abnormalities that may be seen by visual inspection. This 
supports the approach of emphasizing visual inspection for this work. 

Quality control consists of two phases; (i) identification of weld 
discontinuities (some of which may be severe enough to be classified as defects), 
and (ii) assessment or evaluation of the weld based on the observed discontinuities. 
Usually the latter phase results in a pass/fail judgement for the inspected piece. It 
is the conclusion of this work that the first of the above tasks, identification of 
discontinuities, is the most challenging one. It calls for sophisticated image 
processing and image analysis techniques, and frequently ad hoc methods have to be 
developed to identify specific features in the weld image. The difficulty of this task 
is generally not due to limited computing power. In most cases it was found that a 
modest personal computer or workstation could carry out most computations in a 
reasonably short time period. Rather, the algorithms and methods necessary for 
identifying weld discontinuities were in some cases limited. 

The fact that specific techniques were finally developed and successfully 
demonstrated to work illustrates that the general approach taken here appears to be 
promising for commercial development of computerized quality inspection systems. 
Inspection based on these techniques may be used to supplement or substitute more 
elaborate inspection methods, such as x-ray inspections. 
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CHAPTER II 

OBJECTIVES AND SCOPE OF THE WORK 

The main objective of this work was to explore means to incorporate 
computer technology into the process of quality inspection and evaluation, 
particularly for welding applications. The commercial value of such work is of great 
importance, and thus the general approach was oriented towards using technology 
that was affordable and readily available. Visual inspection was selected as the 
primary means for quality assessment, as it has relatively modest equipment 
requirements and a large portion of defects may be indicated by visual inspection. 

Introduction and Objectives 

The role of quality control in the manufacturing processes has become 
increasingly important as national and global competition has intensified in recent 
years. Through technological advances, such as digital computing and signal 
processing, improvements in quality control can be implemented in two key areas. 
First, the process of quality control can be greatly enhanced through the use of 
digital computers for data analysis and records keeping, and through sophisticated 
processing of sensory signals such as acoustic signatures or digitized images. 
Secondly, technological developments can be utilized to automate the quality 
inspection process to a significant degree, and thus make it more efficient and cost 
effective. 

Quality control is an essential part of any manufacturing process as it 
provides the manufacturer and the user of the manufactured item assurance of 
performance and reliability to a specified degree. Expensive failures in large 
projects can frequently be avoided by relatively inexpensive preventive measures in 
quality control. Thus, efforts spent on quality control are frequently well worth their 
cost. 

While most aspects of the manufacturing processes have progressed 
substantially with advances in robotics, electronics, and related fields, the crucial 
task of quality control still relies in most cases on traditional approaches of human 
judgement. Although human quality judgement has served manufacturers quite 
adequately for a long time its drawbacks and limitations have emerged as industrial 
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competition has intensified, complexity of manufactured products and systems has 
increased, and quality requirements and performance specifications have become 
stricter. A number of specific improvements can be accomplished by automating 
the quality evaluation process, as discussed in the following paragraphs. 

Consistency and reliability of the inspection can be improved. It is well 
known that human inspectors can be fairly inconsistent, both individually as a 
person’s judgement criteria varies from time to time, and collectively as one person s 
judgments may vary from those of another one. Although rigid quality inspection 
and manufacturing standards are imposed to ensure adequate product performance 
such schemes will only be as reliable as the people enforcing them. 

Focus on details to an extent unsurpassable by human inspectors is possible 
in automatic systems. While humans usually judge quality by overall appearance 
and at best keep alert against specific details commonly encountered in the 
manufacturing process, a non-human quality controller may be programmed for 
analysis to any specified level of detail. 

Reasoning leading to each judgement and references to previous judgments 
are easily accessible. There will never be doubt to as why a computerized quality 
inspector rejects or accepts an item. All reasoning is performed according to known 
rules or calculations. The information on which a judgement is based is archived so 
it may be accessed at later times by humans or by the computer for its own use and 
reference. 

Speed of evaluation is frequently a concern. By employing appropriate 
hardware and software structures the speed of an automatic quality inspector can 
surpass that of a human inspector. 

Cost, the bottom line, is always significant in manufacturing processes. While 
it can be argued that initially the cost of a computerized quality controller will 
surpass those of human inspectors the cost ratio will be reversed as time proceeds. 

The level of automation in the computerized quality control system may vary, 
depending on the application. In the ideal situation the proposed system would be 
capable of inspection and evaluation in a completely autonomous manner. 
However, a more realistic approach has been taken here in developing a prototype 
system which tackles a specific subset of the entire problem of quality control. 
Furthermore, situations will arise when human interceptions are required to 
evaluate the examined items. The system is therefore designed so as to facilitate 
human interactions. In addition to using the techniques developed for this work for 
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actual quality control tasks, the system may also be used as a training tool for human 
quality inspectors. 

The work carried out during the contract period is to be general enough to 
benefit a variety of quality control applications at NASA, in addition to commercial 
applications, particularly those relating to materials inspection and evaluation. 

Main Results and Suggestions for Furt her Implementations 

The results of this work strongly suggest that visual inspection procedures, 
that so far have mainly been carried out by human experts, may at least to some 
extent be implemented on microcomputers. Quality control has traditionally, and 
for good reasons, been a human-dominated field. Many of the inspection tasks and 
judgments that quality control requires is not easily specified or implemented in 
computer terms, which has left this field to a large degree outside the mainstream of 
computer applications. The rapid development of computer technology in recent 
years, however, has been narrowing this gap between human and automatic 
inspection. 

This report outlines some of the avenues that were explored during the 
contract and the main results. It is shown that a variety of image processing 
techniques exist to enhance and extract certain features from weld images. 
Furthermore, specific methods are illustrated that assign numeric quality indicators 
to welds, based on their visual appearance. 

The work was carried out using readily available computer and imaging 
equipment. The main results consist of two parts, both of which are delivered to 
NASA at the conclusion of the contract. First, the applied research and 
experimentation that was carried out to accomplish the objectives of the work are 
summarized in this report. Secondly, a software package was developed that may 
serve as a template for a commercial visual inspection system. This software 
package is delivered to NASA, as well. 

A proposed continuation of this work will result in a commercial system for 
quality inspection. The software will be modified based on the experience obtained 
with the prototype system, and the hardware will continue to consist of readily 
available components. 
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CHAPTER III 

BACKGROUND - EXISTING TECHNOLOGIES 

A number of techniques are presently widely used for quality inspection of 
metal fabrication and welding in particular. A brief overview of these 
methodologies is given in the following paragraphs. 

Visual Inspection and Other Inspect ion Methods 

While a number of methods are currently used for quality inspections, each 
one has its own area of applicability, advantages, and disadvantages. An extensive 
survey of the state-of-the-art of available methods is vital before the vast task of 
designing and building a comprehensive quality evaluation system is initiated. By 
studying existing testing methods a thorough comparison of the various approaches 
ran be Carried out. The comparison is to illustrate the areas of applicability for each 
of the testing methods, specific strengths, and limitations. This study should enable 
the system designers to program the system to select the optimal testing methods for 
each test situation. 

The testing methods have traditionally been categorized as nondestructive 
tests (NDT) or nondestructive evaluation (NDE) on one hand and destructive tests 
(DT) on the other. As the term indicates NDT methods are those testing methods 
which convey information on discontinuities in the inspected item or material 
without affecting any of its characteristics for usability. This is contrasted to the DT 
methods which require the tested sample to be affected to the extent that it can not 
be put to its intended use after testing. Clearly the NDT methods are preferable to 
DT whenever applicable. This holds particularly for manufacturing of 

small-quantity items and wherever production price is relatively high, due to high 
materials cost, extensive precision, etc. On the other hand DT frequently yields 
information which may not be attainable by NDT. 

Beyond the DT-NDT classification the testing methodologies can be 
categorized according to which defect types they emphasize. Many defects can be 
detected by more than one testing technique and with varying degrees of accuracy. 
Therefore there is a considerable overlap of applicability of the various testing 
methods. By a comprehensive study of the available testing methodologies this 
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research would enable selection of a well chosen set of testing techniques which 
would allow quality inspection for any application in an efficient manner. For 
further clarifications some of the most commonly used quality testing methods will 
be reviewed here. 

Visual Testing (VT) [1] is the most readily applicable testing method for 
human inspectors. Its primary advantage is its relative simplicity as far as 
equipment requirements is concerned. In the simplest applications the inspector 
uses relatively simple tools, such as rulers, micrometers and optical comparators to 
classify and quantify the defects, in addition to a light source and magmfying glasses. 
These tools are, for example sufficient for inspection of many surface defects in 
metals. There are mainly two limitations to this test: First, it only conveys defects 
appearing on the surface of the inspected object, and secondly, it depends heavily on 
the visual acuity of the quality inspector. The types of defects frequently appearing 
on the examined metal or ceramic surface include cracks, nucrofissures and 
sometimes porosity. On welds common defects include overlaps, undercuts, unfilled 
craters, and pock marks. Besides defects directly related to joint itself, visual 
inspection can manifest other deficiencies, such as incorrect joint alignment and 
structure distortion caused by heat, mechanical pressure, or other causes of 
deformation. 

Radiograph Testing (RT) [2] utilizes radiation which penetrates the 
inspected object to allow examination of defects internal to the material. The most 
common radiation types are x-rays and gamma rays. The radiation source is placed 
on one side of the material and a radiation-sensitive detector (a photosensitive film 
or an electronic scanning device) on the other. The photographic approach is the 
most common one, which requires the exposed film to be developed and then 
interpreted by a trained quality inspector. The details of the equipment commonly 
applied for this task and the exact procedures for optimizing image quality will not 
be dwelled upon here. Instead, the applicability and limitation are briefly discussed. 
The radiographic technique identifies most types of material defects, internal as well 
as external. For fusion-joined materials typical defects identifiable by radiographs 
are cracks, porosity, lack of fusion, incomplete penetration, and slag, as well as 
corrosion, fit-up defects, wall thickness and dimensional distortions. The price for 
such generality, however, can be substantial. The equipment needed for 
radiographic inspections consists of radiation sources (x-ray or gamma ray sources 
with appropriate projecting apparatus), films and film holders, film processing 
facilities, and radiation monitoring equipment. Interpreting the radiographs calls 
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for highly trained and skilled personnel. Finally, the radiation is potentially 
hazardous and therefore the exposure requires care and a skilled operator. 

Magnetic Particle Testing (MT) [3] is based on the principle that 
discontinuities in a ferromagnetic piece of metal break up magnetic flux in the metal 
and produce localized magnetic poles at the discontinuities. Magnetic flux is 
conducted through the inspected metal object by placing it in a strong magnetic held 
from an external magnet or by conducting current through the metal. Magnetic 
particles are spread over the surface of the metal and they form patterns along the 
magnetic flux lines. If discontinuities occur at or near the surface they disrupt the 
magnetic flux in the tested object. The disruptions, in turn, reorient the patterns of 
the magnetic particles and thus indicate the defect locations. The equipment 
needed for magnetic particle testing includes magnetizing coils or current prods, a 
current source, magnetic powder, and sometimes ultraviolet light to examine the 
magnetic patterns. While this method is relatively economical, it is limited to 
ferro-magnetic materials and defects located at or near the surface. 

Liquid Penetrant Testing (PT) [4] uses the fact that certain liquids can enter 
voids and cracks and remain there when the surface liquid is removed. Thus, these 
liquid residuals can reveal even the smallest imperfections on the material surface. 
The liquids used for this purpose are either visible in ordinary white light, or in 
ultraviolet light. The testing procedure is simple and requires only the fluorescent 
or dye penetrant liquid, developer and cleaning chemicals and optionally a 
fluorescent light source. The main limitation is that this testing method is only 
applicable to surface voids, such as surface cracks, microfissures, and porosities. 

Eddy Current Testing (ET) utilizes the principle that an electric current will 
flow in any conductor subjected to a changing magnetic field. Application of a 
varying magnetic field to a metal object generates electric currents in the metal, 
whose conductivity usually decreases as the number and size of discontinuities 
increase. An inductor, carrying alternating current, is usually employed for 
generating the magnetic field. Resistance to the eddy currents is reflected into the 
magnetic field and this affects the impedance of the inductor. Thus, the internal 
discontinuities of the metal can be detected by monitoring the inductor impedance. 
Furthermore, the frequency of the inductor current influences the sensitivity of this 
technique to various defect types. The main equipment required for the testing 
includes an induction coil, a current source, magnetic pickup coils, current meters, 
and calibration standards. While the test equipment is usually fairly portable, the 
method is usually limited to surface and shallow discontinuity detection. 
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Ultrasonic Testing (UT) [5] involves directing a high frequency sound beam 
into the test material on a predictable path, which, upon reflecting back from a 
discontinuity in the material, produces a wave that is amplified and usually displayed 
on a graphic device. Thus, the underlying principles are the same as used for radar, 
sonar, and related echo detecting techniques. The ultrasonic waves or pulses are 
generated by electronic means and transmitted to one or more piezoelectric devices 
which transform the electrical energy into mechanical waves. The piezoelectric 
transducers are attached to the tested metal, frequently on a thin layer of liquid, 
which serves as a coupling medium for the acoustic signals. Detection is carried out 
through similar piezoelectric transducers which convert the received sound waves 
into electrical signals. These signals are displayed on graphic devices, such as 
cathode ray tubes. Principally this testing method is very powerful as it can convey 
most defect types and is relatively portable. On the other hand, it requires skilled 
personnel for interpretation of the detected signals. 

Acoustic Emission Testing (AET) [6,7,8] is still considered to be in its early 
stages of industrial use. It consists of the detection of acoustic signals produced by 
plastic deformations or crack formation during loading. Detection is, as in the case 
of ultrasonic testing, carried out by piezoelectric transducers. These signals are 
am plified, passed through signal processing equipment, and interpreted by their 
appearance on a graphic device. AET is useful for continuous surveillance 
inspection, and detection of cracking in loaded structures and during material 
cooling. Fairly sophisticated signal processing and interpretation skills are required 
for this testing technique. 

Of the various testing techniques summarized above, this work has primarily 
been concerned with testing based on visual inspection. As illustrated later in this 
report, images of weld surfaces may be used to evaluate and assess weld quality. 
Image inspection was used as the means for judging the quality of welds from the 
GTAW and GMAW processes. In principle, the same approach may be taken to 
examine VPPA welds. However, the section on VPPAW inspection, later in this 
report, illustrates the use of a laser profile scanner, which may be used to provide 
information on weld surface conditions as well. 

Image Processing and Analysis Techniques 

A summary of image digitizing, analysis and processing techniques will be 
given in the following subsections. This summary should familiarize the reader with 
the approaches discussed relating to this work as well as the terminology involved. 
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The concept of "computer vision" has been separated into three distinct tasks. 
Image digitizing is the first step, where the inspected item is captured by a camera 
and digitized for further digital processing. The second task is preprocessing of the 
image, which generally serves to improve it, eliminate irrelevant or unwanted 
artifacts, and enhance the features that are of importance for the task, i.e., quality 
control. The third step in the computer vision process is analysis of the image, 
which ultimately results in conclusions relating to the features found in the image, 
such as discontinuities or defects in the examined part. 

Image Digitization 

A number of devices exist for capturing images and converting them to 
electronic form for further processing. The light path of most such devices follows 
that of regular cameras, but instead of using a chemical film plane the electronic 
image digitizer uses sensors which convert incident light intensity into electrical 
signals. 

Image sensors based on vacuum tube technology include the image orthicon 
tube, the vidicon tube, the iconoscope tube and the image dissector tube. The 
orthicon tube is a high-quality device in that it offers very stable performance at all 
incident light levels, and because of its qualities it is widely used in commercial 
television. The vidicon is smaller and more rugged than the orthicon, and thus it is 
better suited for portable operations and applications where small size is desirable. 
The iconoscope tube is now mostly of historical interest, and the image dissector 
tube is not widely used either. 

Solid-state sensors are rapidly taking over the functions that vacuum tube 
image sensors used to perform only a few years ago. Broadly, the solid-state image 
sensors can be divided into two subgroups: Charge-coupled devices (CCD’s) and 
charge injection devices (CID’s). Both techniques are based on arrays of 
metal-oxide field effect transistors (MOSFETs), in which incident photons generate 
electrical charges that may be sensed and amplified. In the case of CCD’s each row 
in the two-dimensional array is configured as a shift register with a single output. A 
series of clock pulses is applied to this register and the photon-induced charges are 
read sequentially from the register output. In contrast to this, the CID’s retain the 
electrical charge in its fixed place in the array and the amount of this charge is 
accessed by addressing signals on the semiconductor array. The camera used for 
this work is based on the CCD principle, and the remainder of this discussion will be 
centered on this technology. 
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The resolution of the digitized image depends on the size of the CCD array. 
Usually the array is square (i.e., the number of rows equals the number of columns) 
and typically the number of either rows or columns is an integer power of 2, such as 
256 or 512. Thus, the digitized image, whether obtained through a CID or a CCD, 
may be viewed as a matrix of individual sensor elements, pixels, each of which holds 
a voltage which is related to the incident light intensity. The update rate of the pixel 
values can generally be controlled electronically, and is typically set at 60 Hz (the 

U.S. television system screen refresh rate). 

For digital processing of the image the voltage value of each pixel is digitized 
to a numerical value. This is usually accomplished with electronic hardware 
referred to as image digitizers or frame grabbers. The image digitizer converts each 
pixel in the image to a numerical value which is a function of the brightness or the 
color of the pixel. In the case of color images each pixel may characterized by three 
numbers, each one indicating the intensity of a primary color in that pixel: red, 
green and blue. A white pixel has equal amount of the three primary colors, and its 
brightness is determined by the equal values of these colors. For example, a given 
pixel value may be represented by three color registers in the image digitizer, each 
of which may contain numbers between 0 and 255. With all three register 
containing 255 the resulting pixel color may be bright white, with all three register 
storing the value 127 the color is gray, and with 0 in all register the pixel is black. 
Any color combination can be obtained by mixing the color values in appropriate 
proportions. For example, a green color can be obtained by assigning relatively high 
numbers to the blue and the yellow registers while leaving the red register at a low 
value or zero. An alternative to the red-green-blue (RGB) system is another one in 
which the pixel color attributes are stored as "hue", "saturation" and "intensity" 
(HSI). The intensity component determines the overall brightness of the pixel. The 
saturation is the purity of the pixel color. Highly saturated colors are vivid while low 
saturation results in dull or faded colors. Finally, the hue is the actual color. Using 
either of the above colors schemes, an entire color image can be stored in digital 
memory, where it may be used for further processing. Frequently the image 
digitizer is also capable of converting the digital signals of the stored image to 
analog signals which may be routed to a monitor for viewing. 

Image Preprocessing 

The first step in processing a digitized image is generally referred to as 
preprocessing, which encompasses a number of operations. The purpose of 
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preprocessing is usually to enhance general attributes of the image, such as its 
contrast, average intensity, and thus the preprocessing methods typically operate on 
the entire picture. 

Equalization is an image processing technique, which serves to enhance 
contrasts in the digitized image and neutralize the effects of excessive or inadequate 
lighting. This is useful in real-time welding operation where the light illuminating 
the workpiece is frequently dominated by the arc brightness, which is not readily 
controllable for light intensity. In one relatively simple approach to equalization the 
image is scanned and the maximum and minimum pixel brightness values (intensity, 
in the case of color images) are established. The intensity minimum and maximum 
observed in the image may or may not reach the minimum and maximum that the 
image digitizer offers (e.g., 0 and 255, respectively). Once the original pixel intensity 
range has been determined a new intensity is assigned to each pixel so that the 
entire digitizer range is utilized. Specifically, if and p„ax initial minimum 

and maximum intensities in the scanned image, each intensity, p, will be updated to 
a new intensity value, q, according to the following linear transformation. 

9 * 9niax ( P~PmSn) ^ (Pimx Pmin^ ^ ^ 

where is the maximum pixel intensity level supported by the image digitizer, 
while its minimum intensity level is assumed to be zero. The result of this operation 
is an image that generally has stronger contrasts and is less sensitive to improper 
lighting than the original image. 

Various filtering techniques constitute a host of preprocessing tools. 
Generally, filtering can be considered as a means for transforming the image 
intensities in some way so as to enhance or deemphasize selected features. 

Smoothing or averaging is a relatively simple method for reducing the effects 
of individual pixel abnormalities or "noise" in the image. This operation is 
essentially equivalent to low pass filtering of sequential signals, except in this case 
the signal is the image, which is two-dimensional. With this technique each pixel in 
the image is assigned a new brightness level which is based on its original brightness 
as well as those of the pixels surrounding it. For example, the brightness of each 
pixel in the image may be set as the average of the 3-by-3 array of pixels in which it 
is centered. When this has been repeated for each pixel, the resulting image is 
generally more free of noise 
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A more sophisticated method of equalization than the one described earlier 
in this section is the histogram equalization. The pixels of the image are ordered 
according to brightness values and a histogram of the encountered levels is 
constructed. A bright image will generally reveal large histogram values at high 
brightness values and vice versa for a dark image. Alternatively, a dull image 
contains small histogram values at the two brightness extremes while a high contrast 
image reveals larger numbers at very high and very low brightnesses. Through 
histogram equalization the brightness levels found in the original image are mapped 
so that a histogram of the image after the mapping would result in equal number of 
all brightnesses. The original image may produce a histogram h(p) and the 
his togr am of the equalized image may be denoted ss g((j), where p and <j represent 
brightness levels. Clearly, the number of any given set of original pixels with 
brighmess in the range p to p+dp equals the number pixels with corresponding 
brightnesses between q and q+dq after the histogram mapping. In other words: 

h(p) dp = g(q) dq (2) 

Furthermore, if the new histogram, g(q), is to be uniform, then: 

g(q) dq = ff/M (3) 

where N is the number of pbcels in the image and M is the number of brightness 
levels. The histogram equalization transform is therefore found as: 

q(p) = (M/N^)2 h(s) ds (4) 

0 


Image Analysis and Feature Extraction 
Template matching is a relatively simple operation that facilitates the 
detection of a specified feature in the image. The feature that is sought must be 
well defined and known in terms of dimensions and brightness levels. Usually 
template matching is accomplished through correlation. Basically, the correlation 
procedure is carried out as follows. The sought feature, e.g., a specified item in an 
image, is defined on a grid of pixels, a template, that may be moved around the 
analyzed image. For each location of the template in the image the difference 



13 


between each template pixel brightness and the corresponding image pixel is 
calculated, squared, and the sum of all such squares over the entire template is 
determined. This sum is an indicator of the overall difference between the template 
and the part of the image located under the template. Specifically, the sum is 
expressed as: 

D(xJ = 2 [f(x)-t(x-z)P (5) 

X 

Here, f(x) is the image, which consists of pixels x = (Xi,Xi), and t(x-^ is the 
template, displaced to location y = (YisYz). The sum D(y) varies with y and is 
taken over all of the template pixels. The squared term may be expanded as: 

[f(x)-t(x-YjP = f^(x) - 2f{^t(x-Y) + t^(x-Y) (6) 

Clearly, is independent of y, and thus it is constant. Similarly, the template, 
tHx-Y). is by definition fixed, which leaves the center term as the variable 
component. The variable component is proportional to the cross correlation 
between f{x) and t(x-Y), which is defined as: 

Rn(Y) = 2 [f()0 t(x-Y)l (V 

X 

and is maximized when the image under t() is identical to t() itself. Thus, 
correlation can be used to search for a known item or feature in the image and 
determine its location as well. 

Edge detection is an operation that is concerned with locating areas of 
abrupt brightness variations in the image, such as at edges or boundaries between 
objects. A simple edge detection operator is one that replaces each pixel with the 
magnitude of the brightness gradient at that pbcel. For example, each pixel 
brightness p(x), for a pixel located at x=(xl,x2), may be replaced by the quantity 
G(x), determined as: 

G^(x) = [p(Xi+l,X 2 )-p(X^,X 2 )P + [p(Xi,X 2 +l)-p(Xi,X 2 )P (8) 

Edge detection is useful for enhancing artifacts in the image, such as cracks 
or other sharp features. Edge detection needs to be used carefully, however, as it 
tends to highlight arbitrary noise pixels in the image, as well. 
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CHAPTER IV 

IMAGE ANALYSIS AND EVALUATION FOR QUALITY CONTROL 

While the previous chapter outlined some of the general aspect of image 
processing and common image processing techniques, this chapter will be devoted 
to specific techniques relevant to image processing for quality control. 

As stated earlier, evaluation of metallic welds has been selected as a specific 
quality control application for this work, although other types of materials 
inspection may benefit from the work as well. Evaluation of weld quality takes a 
number of factors into account and is, when taken to its ultimate degree, a fairly 
sophisticated process. In this work a number of specific features of the weld bead 
were selected to represent the overall weld quality. The set of examined features 
may not always contain every possible defect type, but some of the most common 
characteristics separating "good" welds from the "bad" ones were exarmned for this 
purpose. 


Feature Enhancement 

Most visual weld defects are characterized as sharp artifacts in the image. 
This includes cracks, weld bead ripple patterns, spatter and porosities. Edge 
detectors were therefore important tools for enhancing these weld features. 

The gradient operator, described earlier, was used for most edge detection 
tasks. Figures 1 through 3 illustrate examples of edge detection applications. 
Figure 1 shows a metallic piece with a crack propagating from its edge. In figure 2 
the crack has been enhanced by a magnetic particle test. Applying the edge 
detection operation on this image highlights all sharp features, resulting in a much 
enhanced view of the crack, as shown in figure 3. 



Figure 1. An unprocessed image showing a crack emerging from a machine 


part. 




Fimre 2. An image of the CTack in the workpiece after magnetic particle 
material has been applied. 



Figure 3. The image of the aack after edge detection algorithm has been 
applied. 


Analysis Using Artificial Neural Networks 
Artificial neural networks (ANN’s) have gained prominence recently among 
researchers of various practical problems which have not been solved adequately by 
traditional methods. As the name implies, these networks are computer models of 
the processes and mechanisms that constitute biological nerve systems, to the extent 
that they are understood by researchers. 

Interest in artificial neural networks has grown rapidly in the past years as 
they have been envisioned by many as an alternative to artificial intelligence (AI) as 
a way to attack problems that "traditional" AI methods have been unable to solve. 
Most of the growing interest in neural networks has surfaced since the mid-1980’s 
and by now, at the beginning of the 1990’s, neural network conferences are arranged 
by professionals from such diverse areas as science, banking and financing. 
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engineering, snd the military. The widespread interest in neural networks stems 
from their ability to solve complex problems which so far have proved difficult for 
competing methods. 

In addition to the mere ability to conquer complex problems, neural 
networks are attractive in view of their high execution speed and their relatively 
modest computer hardware requirements. Beyond the variety of personal computer 
an d workstation implementations that are available commercially, several efforts 
are presently underway to take advantage of the inherent parallelism of neural 
networks. The calculations taking place within the networks are not generally 
carried out in a serial fashion as they can be conducted in parallel. Although 
present computers implement neural networks very efficiently, based on the 
traditional, serial, Von-Neuman architeemre, emerging parallel computers will be 
ideally suited for neural network implementation. 

Neural networks were used in this work to analyze weld profiles and 
determine the location of certain features, such as the bead center, undercuts, the 
edges of the bead, etc. This technique has been demonstrated to be particularly 
effective in analyzing the VPPA weld profiles. Specific examples are given in 
chapter V, and further discussion on artificial neural network is provided in 
appendix A of this report. 


Rule-Based Expert Systems 

The use of rule-based expert systems for weld quality assessment was 
examined during the course of this work. A number of approaches exist to process 
the information obtained from the image processing routines and expert systems 
were studied for this purpose. Among the questions that arise are what constitutes 
an "acceptable" workpiece or specimen, and to what extent defects can be tolerated. 

A number of software tools were surveyed for the purpose of making 
PASS/FAIL decisions. A relatively logical candidate was the CLIPS (C Language 
Integrated Product System), which has been developed by NASA at the Johnson 
Space Center [9], and is available to NASA contractors for work such as this one. 
This system allows the user to specify a collections of general "rules" and define a 
collection of "facts" relating to the problems at hand, as shown in figure 4. The 
CUPS "inference engine" is used to iteratively compare the rules against the facts, 
which in turn derives sets of new facts, and thus the system continues until no 
further deductions can be made. 



Figure 4. A schematic diagram of the CLIPS rule-based expert system. 

A small expert system was constructed to examine the feasibility of 
determining weld quality based on observed defects. Representative samples of the 
facts list and the rule base are listed below. 

Facts ; 

(deffacts qc-facts 

(defect-types longitudinal -crack crater-crack transverse-crack)) 
Rules: 


Rule 1 
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(defrul e 1 ongi tud i nal +crater-crack 

(longitudinal -crack) & (IF....) 

(crater-crack) 

=> 

(assert (quality =4))) (THEN...) 


Rule 2 

(defrule longitudinal -crack 
(longitudinal -crack) 

(assert (quality - 6))) 


Rule 3 

(defrule crater-crack 
(crater-crack) 

=> 

(assert (quality =8))) 


Rule 4 

(defrule transverse-crack 
(transverse-crack) 

=> 

(assert (quality = 5))) 


Rule 5 

(defrule porosity 
(porosity) 

=> 

(assert (quality =6))) 


The above system and other similar expert systems were evaluated and 
generally worked as expected. It was found, however, that the main problem 
associated with the weld evaluation was not in the reasoning process (which the 
expert system is mostly concerned with) but rather in the image analysis process. 
The evaluation schemes that were finally selected and demonstrated to yield 
acceptable results are discussed in chapter V. 
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CHAPTER V 
QUALITY EVALUATION 

The task of weld quality evaluation is not trivial, even for the experienced 
quality inspector. This is particularly true when it comes to specifying in 
quantitative terms what attributes of the weld affect its quality and to what extent. 
Different types of discontinuities have been categorized for this purpose, such as 
cracks, porosity, undercuts, microfissures, etc. The task of quantifying the presence 
of these discontinuities and relate this to the overall quality of the weld is not easily 
specified, and the task is even more difficult when it comes to coding them in terms 
of computer operations. 

An overview of the quality evaluation process is outlined in figure 5. The 
block diagram contains the functions which are carried out in the quality evaluation 
process. Examples of particular types of each function are shown on the left. All or 
none of these at each step may be used in any particular evaluation. The workpiece 
to be inspected may optionally be treated by a variety of standard pre-inspection 
methods, such as magnetic particle treatments or dye penetrant application. The 
piece is appropriately illuminated with white light, colored light, structured light, or 
laser illumination, and an image is captured and digitized. Various image 
enhancement procedures are available to extract features from the image, such as 
enhance cracks, porosities or other artifacts. The image processing techniques 
applied to the digitized picture include integration, thresholding and other standard 
image processing methods. Finally, an evaluation of the inspected weld is carried 
out, resulting in a quality assessment of the inspected workpiece. 
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Figure 5. An overview of the quality evaluation process. 

The approach taken for this work was to demonstrate the use of computer 
analysis to evaluate welds based on global characteristics of their appearance. 
Indicators of weld quality include the general regularity of the weld bead, the 
relative spacing and regularity of the solidified pool ripples along the bead, and the 
symmetry of the weld. Other approaches may be taken to the issue of quality 
control for welding, but the ones used proved to produce consistent and reliable 
indicators of weld quality when compared against objective human quality 
judgments. 

Due to the different characteristics of various welding processes and the 
difference in weld appearance, it was found necessary to apply different evaluation 
strategy to each process. This was done for three common welding processes. The 
Gas Metal Arc Welding (GMAW) process is a very common welding process and 
widely used in various commercial welding applications. One of its advantages is 
the high metal deposition rate, which makes it attractive for high-quantity 
applications. A continuously-fed wire serves both as the electrode maintaining the 
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arc and as a filler metal which melts into the welded joint. Inert gas flow surrounds 
the arc and the molten pool. The Gas Tungsten Arc Welding (GTAW) process is 
widely used for precision work. Unlike the GMAW process, GTAW uses a 
non-consumable tungsten electrode and filler wire is fed separately into the joint. 
Finally, the Variable Polarity Plasma Arc Welding (VPPAW) process is of some 
interest here, particularly because it is increasingly used by NASA and commercial 
users. It is similar to the GTAW process with two notable differences. First, two 
sources of inert gas are used; one serves to form the plasma arc, while the other 
serves primarily for shielding the arc and the pool from the atmosphere. Secondly, 
the arc current polarity is periodically switched between the nominal 
electrode-negative mode and brief electrode-positive intervals which serve to 
cleanse the oxide layer that forms on metals such as aluminum. 

Because of the substantial differences between these welding processes, the 
evaluation strategies will be outlined separately in the following sections. 

Quality Evaluation for the Gas Tungsten Arc Weldi ng Process 

A good GTAW weld is typically characterized by a regular weld bead 
surface, as illustrated in figure 6. The surface is free from any indications of 
cracking, porosities, undercuts or other discontinuities. Furthermore, a regular pool 
ripple pattern forms along the weld bead. This pattern tends to get irregular and 
distorted if the welding speed is out of the range that is otherwise suitable for a weld 
pool of the given size. For comparison, figure 7 shows a typical low quality weld, 
which would hardly pass any reasonable inspection process. 
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Figure 6. A photograph of a good quality GTA weld bead. 



Figure 7. A photograph of a low quality GTA weld bead. 
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As discussed earlier, the ripple pattern serves an important role in the 
approach taken here to weld quality evaluation. Two primary criteria were selected 
for evaluation of the ripple pattern. First, increased regularity generally indicates 
increasing quality. Skipped ripples or clustering of ripple formation are good 
indicator of poor welding. A second, and somewhat independent, measure of ripple 
pattern quality is the average spacing between the ripples. Under certain 
circumstances relatively regular wide spacing may be apparent, caused by too high 
travel speed and thus "regular" skips in the ripple pattern. In a summary, a high 
quality GTA welds may be characterized by (i) regular ripple spacing and 
(ii) appropriate ripple/bead width ratio. 

The ripple pattern can be analyzed through the "integrate" menu selection in 
the image processing software developed for this work. This operation consists 
essentially of two tasks. First, the bead width has to be determined from the 
inspected image, and then the spacing of the ripples along the weld bead must be 
determined. Generally, the evaluation technique normalizes the ripple spacing to 
the bead width. By doing this, the focal length of the camera, and the working 
distance between the camera and the workpiece do not have to considered in 
analyzing the ripple spacing. The bead width may be estimated with the horizontal 
integration procedure of the image analysis software. This procedure is graphically 
illustrated in figure 8. The horizontal dimension of the sampled image measures 
pixels and the vertical dimension is Yf^-l pixels. For the GTAW process the 
window of concern extends between the bottom and the top of the image, and its 
width may be varied by the user. Integration is carried out along the horizontal 
(x-axis) direction of the image. The result is a pattern that represents the weld 
profile as seen along the weld axis, as shown in the figure. 
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Figure 8. The horizontal integration procedure illustrated. This technique 
is used in analysis of GMA weld quality. 

Actual image data obtained through the integration process is shown in 
figure 9. The shape of the weld profile is clearly recognized, with a surface peak at 
the center and the bead edges on both sides. Although the general pattern is clearly 
recognized by the human eye, a computer recognition of the bead features is not 
straightforward. Using filtering, however, the signal can be improved. Figure 10 
shows the same pattern, filtered by 21-pixel runmng averaging, which presents a 
better signal for computer processing. The bead profile is used here to estimate the 
bead width to establish a reference for the bead ripple spacing. 
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Figure 9. Image scanning data representing a GTA weld profile. 
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Figure 10. Image scanning data representing a GTA weld profile, filtered to 
aid further computer analysis. 
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A 16-line high window was used for the horizontal integration procedure, and 
the brightness intensities of the 16 pixels of each horizontal position were integrated 
to obtained a somewhat smoothed image of the ripple pattern along the weld bead. 
The vertical integration technique is graphically illustrated in figure 11. The 
horizontal dimension of the sampled image measures Xf^-l pixels and the vertical 
dimension is K„-J pixels. For the GTAW process the window of integration is 
generally placed along the centerline of the image. The ripple pattern is obtained 
from the specified window, as shown at the bottom of the figure. 
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Figure 11. The vertical integration procedure illustrated. This technique is 
used in analysis of GTA weld quality. 

A variety of methods were investigated to yield information regarding the 
weld ripple spacing characteristics. The Fourier transform is a natural candidate for 
this task. A regular ripple pattern with consistent spacing should yield a Fourier 
spectrum with a single, distinct peak at the frequency corresponding to the ripple 
spacing. Furthermore, a narrow peak should indicate small variations in ripple 
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spacing, while a broader peak would result from larger variations in the pattern 
spacing. The Fourier transform was taken of several ripple pattern and examined. 
Typically, the Fourier spectra turned out to be very irregular, as figure 12 illustrates. 
This spectrum was obtained from a weld which had been evaluated to be of a 
relatively good quality. Despite a dominating spectrum peak, corresponding to the 
dominating ripple frequency, the irregularities are severe enough to produce 
additional peaks of significant amplitude. Typically, the power spectrum was much 
less clear for welds of less quality, and frequently no spectrum peak could be 
determined to be the dominant one. 



Various smoothing techniques were tried to clarify the Fourier spectrum, in 
order to enhance further the primary frequency peak. None were found to work 
satisfactorily over a considerable range of weld qualities. In addition to Fourier 
analysis, autocorrelation techniques and other related methods were investigated, 
but none yielded reliable, consistent results. 
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Following the investigations outlined above, a new technique was adopted 
for weld evaluation. The selected approach to computerized GTA weld quality 
evaluation is summarized in figure 13. The video image of the bead surface is 
sampled, digitized, and stored in computer memory. Two aspects of the weld bead 
are quantified: the bead width and the relative density (spacing/bead width, or 
frequency) of the ripple pattern along the weld. These quantities are obtained by 
analysis of windows from the image that have been integrated along the horizontal 
and vertical dimensions, respectively. The ripple pattern frequency and the variance 
of this frequency are used to assess the overall quality of the GTA weld bead. 



Figure 13. An overview of the process developed for evaluating GTA weld 
quality. 

To demonstrate further the technique used for quality evaluation, four GTA 
welds were arbitrarily selected and evaluated on a scale from 0 to 10 by a human 
expert. The GTAW quality evaluation process was based upon observation of three 
major criteria. The first was regularity of spacing of the ripple pattern on the 
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surface of the weld. Uniformity in this aspect of the weld’s appearance was 
considered essential to assignment of a high quality rating. However, regularity 
alone was not sufficient to achieve a high quality rating. The relative frequency of 
the ripple spacing compared to the bead width was also a significant factor in the 
assessment process. A low relative frequency (approximately 1/1) was given a lower 
rating than a high relative frequency (approximately 10/1). The third factor in the 
rating process was the presence or absence of faults such as cracks, porosity, abrupt 
changes in lateral or vertical cross section, or other observable defects. 

For the GTAW test samples, a graphically produced "ideal" weld, with 
perfectly uniform ripple spacing, with excellent relative frequency, and with 
perfectly uniform geometry, was given the score of ten. A real weld, with excellent 
uniformity, excellent relative frequency, and no significant geometrical anomalies, 
was given a rating of nine. This weld was shown in figure 6. A weld with good 
uniformity, but with very poor relative frequency, and severe variations in geometry, 
was given the rating of two. This weld was shown in figure 7. Between these 
extremes, other welds were inspected and given subjective ratings based upon the 
presence or absence of the factors noted above. The welds will be examined further 
as the general technique is illustrated. 

Once the ripple patterns had been digitized they were filtered for noise 
reduction. The filtering method was a 5-pbcel running average. In other words, the 
intensity of pixel k was replaced by the average of the intensities from pixels k-2 
through k+2. The original ripple patterns and the filtered versions are shown in 
figures 14 through 21. The human-assigned weld quality is indicated in the name of 
the file storing each ripple pattern. For example, file GTAW2.SCN stores scaimed 
ripple data from a GTA weld that was evaluated of quality 2. A number of welds 
were originally examined for this purpose and whenever more than one weld was 
assigned the same quality score by the human expert, successive welds were 
designated by letter suffixes, such as 2A, 2B, etc. Thus, the welds illustrated here, 
GTAW2, GTAW4B, GTAW6 and GTAW9B, represent a broad range of weld 
quality levels. Weld quality level 10 was simulated by use of a graphically produced 
weld image with perfect uniformity of width and ripple spacing. 
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Figure 14. Scanned ripple pattern of weld GTAW2 (human quality 
evaluation = 2). 
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Figure 15. Filtered ripple pattern of weld GTAW2 (human quality 
evaluation = 2). 








33 



Figure 16. Scanned ripple pattern of weld GTAW4B (human quality 
evaluation = 4). 



Figure 17. Filtered ripple pattern of weld GTAW4B (human quality 
evaluation = 4). 
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Figure 18. Scanned ripple pattern of weld GTAW6 (human quality 
evaluation = 6). 
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Figure 19. Filtered ripple pattern of weld GTAW6 (human quality 
evaluation = 6). 
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Figure 20. Scanned ripple pattern of weld GTAW9B (human quality 
ition = 9). 


Figure 21. Filtered ripple pattern of weld GTAW9B (human quality 
evaluation = 9). 
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Once the ripple patterns had been obtained through the image digitization 
system and filtered accordingly, a simple algorithm was used to determine the 
pattern frequency and regularity. First, the spacing between adjacent ripples could 
be determined by measuring the distance between adjacent ripple peaks in the 
filtered graph. A pixel was determined to be at an intensity peak if both the pixel 
immediately preceding it and the pixel immediately following it had lower intensity 
values than the examined pixel. Using this classification, the ripple spacings were 
determined. To eliminate the effects of differing weld sizes, the ripple spacing 
values, originally determined as number of pixels, were normalized to the width of 
the weld (in number of pixels). Referring to the two quality criteria stated earlier, a 
good weld is indicated by relatively dense ripple pattern (small spacing/bead width 
ratio) and pattern regularity (small standard deviation of the measured ripple 
spacings). Using these two criteria, the rule used here for quality evaluation was: 

Q = 1 / [(average relative spacing) 

+( standard deviation of spacing)] (9) 

where (? is a positive numerical indicator of weld quality and increases as the quality 
of the examined weld improves. The average spacing between adjacent weld ripples 
is determined as described above and the standard deviation is found similarly. 

The quality values of the examined welds were calculated and graphed 
against the quality value assigned by the human inspector. This is shown graphically 
in figure 22, where the human-assigned quality values are assigned to the horizontal 
axis and the vertical axis represents the calculated quality indicators. It should be 
noted that exact numerical correspondence between the human indicator and the 
calculated one is not sought. Of main concern is that the graph is monotonic, which 
ensures that the human and the computer agree on comparison between individual 
welds. In other words, if the human evaluates weld B to be of better quality than 
weld A, the computing algorithm should agree with this judgement. 
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Figure 22. Comparison of the calculated quality indicator and the 
corresponding human quality evaluations. 

As seen in figure 22 the algorithm developed for evaluating GTA welds is 
consistent and generally agrees with the quality evaluations of the human expert. 
The technique developed for this work resulted from the desire to develop 
algorithms which strongly parallel the human reasoning process in evaluating the 
ripple patterns on the weld surface. Interestingly, the technique selected for this 
evaluation is relatively simple in that it does not rely on overly complicated 
mathematical or signal processing techniques, such as spectral analysis, 
autocorrelation or cepstrum analysis, all of which were used in preliminary 
investigations for this study. Yet its results surpass the preliminary findings 
suggested by those advanced methods. As a conclusion, this technique is an 
illustrating testimony to the value of exploring the human evaluation process and 
using it directly to supplement or replace more orthodox signal processing 
techniques. 
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Quality Evaluation for the Gas Metal Arc Welding Process 

A good GMA weld is not as strongly characterized by a dense, regular ripple 
pattern, as the GTA weld is. In many cases a perfectly good GMA weld does not 
show a distinct trace of surface ripples. Thus, an alternative strategy was developed 
to evaluate welds from the GMAW process. 

Figures 23 and 24 show good quality and poor quality GMA welds, 
respectively. Generally, good quality GMA welds are uniform and contain little or 
no artifacts on the bead surface. Furthermore, the bead width is relatively uniform 
along the length of the bead. In contrast to this, a poor GMA weld typically shows 
signs of spatter, irregular fusion line with the parent metal on the sides of the weld, 
and various discontinuities on the weld surface. The human quality assessment was 
based on the tendency of a good GMA weld to have a very uniform geometry and 
surface appearance along the longitudinal axis of the weld. A weld which is 
excellent in these characteristics was given the rating of nine. This weld is 
illustrated in figure 23. A weld which was performed in the presence of severe 
surface contamination, which resulted in significant variations in the height, width, 
and contour of the weld bead was given a rating of two. This weld is shown in 
figure 24. A GMA weld performed without shield gas, which resulted in severe 
geometry variation, along with extensive porosity, was given a quality rating of one. 
Between these extreme cases of quality, a series of GMAW welds were evaluated 
and assessed a quality rating according to the degree of bead uniformity (geometry 
and appearance) along the weld axis. 
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Figure 23. A photograph of a good quality GMA weld bead. 



Figure 24. A photograph of a poor quality GMA weld bead. 
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The general approach used to evaluate GMA welds was to examine the 
variance of weld profiles, sampled across the weld bead in various locations. A high 
quality weld would generally yield small variance in the weld profiles, while low 
quality weld profiles would vary substantially, as irregularities and various 
discontinuities are encountered in the distinct profile scans. Unlike the GTAW 
analysis, the GMAW evaluation is primarily based on the horizontal integration 
technique, which was discussed earlier, in the section on GTAW evaluation. The 
computerized GMA weld quality evaluation technique is summarized in figure 25. 
The video image of the bead surface is sampled, digitized, and stored in computer 
memory. Ehiring system development, two aspects of the weld bead were analyzed. 
The main aspect was the variance of the weld profile, as examined in different 
locations along the weld. A segmented horizontal integration technique was 
implemented for this. The second aspect of the weld inspection was analysis of 
brightness variations along the bead edges. Generally, it may be anticipated that a 
good weld has more uniform edges than a bad one. For this purpose a full 
horizontal integration was used to determine the edges of the bead, and then a 
vertical integration was used to examine the regularity of the bead edges. This 
second aspect of the weld inspection did not always turn out to be as reliable as the 
first one, i.e., the segmented horizontal integration followed by the variance analysis. 
Nevertheless, it is mentioned here to illustrate the various approaches tested for 
evaluation purposes. 
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VIDEO tIONAL 
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(QMAW PROCEDURE) 


Figure 25. An overview of the process developed for evaluating GMA weld 
quality. 

As for the GTAW evaluation, the general approach for the GMAW process 
is illustrated through examples. In this case five welds were examined, each one 
assigned a quality number in the range from 0 to 10 by a human inspector. The 
selected welds were evaluated with quality indicator values of 2, 4, 6, 8 and 9, and 
these numbers are embedded in the file names containing image data from these 
welds. For example, image file GMAW6H1 contains data from a GMA weld 
evaluated by a human to be of quality 6. The file contains a histogram that was 
scanned horizontally by the image analysis system (thus the letter H), and the file 
represents scan number 1 of this weld. The data from this file is illustrated in 
figure 26. As explained in the chapter on the software system, horizontal integration 
essentially averages to a specified degree the pixel intensities in the image. The 
figure therefore illustrates a filtered image of the weld profile, as viewed along the 
weld axis. 
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File GMA6H1 ("Quality" = 6 ) 
Scan No. 1 Across Profile 



Figure 26. Horizontal integration scan No. 1 of weld GMA6. 

While figure 26 shows one scan across the bead, the profile variance along 
the bead length was assessed by taking a total of 5 scans, each of 16 pixels width. 
Figures 27 through 31 show the result of this where, in each case, the same weld has 
been scanned at five different locations and the resulting scans are all graphed 
together in the same illustration. The similarity between the individual bead 
profiles obtained from the same bead, at different locations, may be used as an 
indicator of the uniformity of the weld and thus its overall quality. 
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Files GMA2H1-GMA2H5 ’* = 

Scanned at 5 Locations Across Profile 



Figure 27. Horizontal integration scans of weld GMA2. 


Files GMA4H1-GMA4H5 ("Quality’' = 4) 
Scanned at 5 Locations Across Profile 



Figure 28. Horizontal integration scans of weld GMA4. 




Relative Brightness Relative Brightness 
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Figure 29. Horizontal integration scans of weld GMA6. 


Files GMA8H1 -GMA8H5 ("Quality" = 8) 
Scanned at 5 Locations Across Profile 



Pixels, Across Bead 


Figure 30. Horizontal integration scans of weld GMA8. 
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Files GMA9H1 -GMA9H5 ("Quality" = 9) 
Scanned at 5 Locations Across Profile 



100 150 

Pixels, Across Bead 


250 


Figure 31. Horizontal integration scans of weld GMA9. 

At least two approaches may be taken to quantify the similarity (or lack 
thereof) of the five scans. One is to measure the envelope of the ensemble of scans 
at each location, i.e., measure the difference between the maximum and the 
minimum intensities. The other approach is to calculate the root-mean-square 
(RMS) variance of the five points of each location. Both methodologies were 
implemented and tested. Generally, it was found that they yielded similar results. 
Whenever the envelope measure was large, so was the RMS measure, and vice 
versa. The calculated envelope differences and RMS variances were determined for 
each of the tested welds and the results are illustrated in figures 32 through 41. 




RMS Brlj^ness Spread, Relative ^ (Max-MIn) Brightness, Relative 
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Files GMA2H1 -GMA2H5: 
Envelope Spread, Avg=0.328 



Figure 32. The envelope difference of scans across weld GMA2 (human 
y evaluation = 2). 


Files GMA2H1 -GMA2H5: 
RMS Spread, Avg=0.121 



Figure 33. Root-mean-square variance of scans across weld GMA2 (human 
quality evaluation = 2). 





Files GMA4H1 — GMA4-H5: 
Envelope Spread, Avg=0.146 



Figure 34. The envelope difference of scans across weld GMA4 (human 
quality evaluation = 4). 



Figure 35. Root-mean-square variance of scans across weld GMA4 (human 
quality evaluation = 4). 



RMS Brightness Spreod, Relative (Max-Mln) Brightness. Relative 
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Files GMA6H1 -GMA6H5: 
Envelope Spread, Avg=0.115 



PIxals, Across Bead 


Figure 36. The envelope difference of scans across weld GMA6 (human 
y evaluation = 6). 


Files GMA6H1 -GMA6H5: 
RMS Spread, Avg=0.046 



Figure 37. Root-mean-square variance of scans across weld GMA6 (human 
quality evaluation = 6). 





RMS Brightness Spreod, Retative (Max-MTn) Brightness. Relative 
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Flies GMA8H1 -GMA8H5: 
Envelope Spread, Avg=0.110 



1 

Figure 38. The envelope difference of scans across weld GMA8 (human 
y evaluation = 8). 

Files GMA8H1 -GMA8H5: 
RMS Spread, Avg=0.04-4 
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Figure 39. Root-mean-square variance of scans across weld GMA8 (human 
quality evaluation = 8). 
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Figure 40. The envelope difference of scans across weld GMA9 (human 
quality evaluation = 9). 



Figure 41. Root-mean-square variance of scans across weld GMA9 (human 
quality evaluation = 9). 
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The average of these differences (envelope or RMS measures) yields a 
indicator of the overall dissimilarity of scans within each weld. The average of each 
measure is printed at the top of each of the preceding graph. It may be seen that 
the envelope difference and the RMS variance generally decrease for improved 
weld quality, as expected. Rather than using only one of the variance measures, 
both may be combined into a quality indicator of the form: 

* J / [(avg envelope spread)+(RMS variance)] (10) 

This quality indicator is illustrated in figure 42 as a function of the 
human-assigned quality values. As for the GTAW case, match of the actual values 
is not important, but rather consistency so that comparison between any two welds 
yields the same result as a human would come to. 
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Figure 42. The GMAW computer quality indicator as a function of the 
human-assigned quality values. 
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Duality Evaluation for the Va r iable Polarity Plasma Arc Wglding Prpcgss 

The third welding process of interest for this work is the VPPAW process. 
This is a relatively new welding process, which is increasingly used by NASA as well 
as in conunercial applications. 

The VPPA weld appearance is to a large extent similar to that of the GMA 
weld. Unlike the GTA weld, which is characterized by a ripple pattern along the 
bead, the VPPA bead surface usually has negligible or no surface ripples at all. The 
weld edges, however, may become irregular under certain conditions, a 
characteristic that is shared by the GMAW process. Actual VPPA weld samples 
were not available for this work to test the GMAW quality assessment method on 
VPPA welds, but the method may be assumed to be equally applicable to both 
processes. An alternative approach is presented here, however, which is based on 
earlier research [10] in which the VPPA weld bead profile was analyzed with a laser 
scanner for symmetry and other features that may be used to assess weld quality. 

Figure 43 illustrates the physical components of the welding process. The 
welding torch moves along the joint at a constant speed. In this case a laser scanner 
is used to obtain the weld profile information. Unlike a camera, which essentially 
records the light patterns resulting from reflections and shadows from the bead 
surface, the laser scanner determines the distance to the bead surface. Thus, each 
scan of the laser yields the true height profile of the bead along the scan across the 
weld. Through repeated scanning across the bead, as the scanner moves along the 
weld with the torch, this method may be used to record the true three-dimensional 
topography of the bead. For the quality control technique illustrated here, however, 
individual scans are obtained and processed separately. Each profile scan is routed 
to a neural network that analyzes the profile and detemunes specified 
characteristics, such as its symmetry, dimensions, and other desired features. These 
characteristics are used finally to assess the quality of the weld. 
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Figure 43, The VPPAW process and the components used to scan and 
analyze the weld profile. 

The weld quality obtained with the VPPAW process may be assessed with a 
similar method as the one used for GMAW, where scans across the weld bead were 
compared for irregularities. Figure 44 shows a photograph of a typical high quality 
VPPA weld cross section. Sequential laser scans of the top of this weld can be used 
to reconstruct the three-dimensional weld surface profile, resulting in similar set of 
data as used for the GMAW process. 










Figure 44. A photograph of the cross section of a ^ical VPPA weld. 

The general VPPAW evaluation process is outlined in figure 45, assuming 
tha t the laser scanning method is used. The profile data fi’om a single laser scan is 
sampled and recorded. A neural network is trained off-line with a part of the 
available data, so that it may analyze general features of the bead profiles. The 
results yielded by the neural network are used to evaluate the VPPA weld for 
various quality measures, such as bead symmetry, dimensional regularity, etc. Any 
number of quality assessment criteria may be constructed to utilize the results from 
the neural network. For example, the offset of the crown peak fi'om the true center 
between the undercuts or the bead edges may be used as an indicator of degraded 
quality. The same may be applied to mismatch in the depths of the undercuts. 
These features will be discussed further in the following paragraphs. 



55 



Figure 45. The procedure for VPPA weld evaluation. 

Using this approach, several scans may be made across the VPPAW profile 
and the resulting profile data may be analyzed for variance, which in turn may be 
used as an indicator of weld irregularity. In addition to this technique, which was 
discussed in the section on GMAW weld evaluation artificial neural networks may 
be used to determine the locations of key features of the weld profile, such as its 
crown, the undercuts and the bead edges. The neural network system used to 
analyze the digitized weld profiles is outlined in figure 46. The height values of the 
digitized profiles were entered into the 60 inputs of the network, which was trained 
to locate the five features of the profile and present these locations at the outputs. 
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Figure 47 illustrates a typical cross sectional weld profile. The objective of 
the neural network was to find the locations of the weld crown, the left and right 
undercut, and the left and right weld boundaries. To the human observer, these 
features are relatively clear, with the possible exception of the weld boundaries. 
The locations of these features are indicated on the figure with the arrows 
marked "N”. For comparison, an alternative system, presently used at NASA, 
located the same features as shown with the arrows pointing from below the profile 
trace. In this case the present system failed in correctly identifying the right weld 
boundary, while the neural network succeeded. In most cases the present system 
performed satisfactorily, while occasionally it gave incorrect results. The errors 
usually occurred in determining the weld boundaries, and these errors could be 
significant. For comparison, the neural network performed somewhat better, on the 
average, and the worst case errors were always much smaller than those of the 
present system. 
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Figure 47. A scanned weld profile, with the weld crown, undercuts, and 
parent metal boundaries indicated by the neural network and by the traditional 
system. 

A set of 60 bead profiles were used to test and train the neural network. Of 
the total of 60 profiles, 30 were used for training and 30 were used for testing the 
trained net. 6 of the testing profiles were sampled for error analysis, and the results 
are listed in table 1. As shown there, the errors in determining the locations of the 
bead crown, undercuts and parent metal boundaries are typically less than 1%. 




58 


TABLE 1 

Sample results obtained with the weld profile anal)^sis neural network. These six 
weld profiles were not used in original training of the network. 



L-Par 

L-Ucut 

Crown 

R-Ucut 

R-Par 

RMS 

Profile No.: 

Err: 

Err: 

Err: 

Err: 

Err: 

Err: 

134.17 

-0.5% 

0.6% 

-0.4% 

-0.2% 

1.8% 

0.9% 

148.17 

0.6% 

0.7% 

-03% 

-0.1% 

-1.1% 

0.7% 

302.17 

0.2% 

0.6% 

1.0% 

0.8% 

0.0% 

0.6% 

316.17 

0.4% 

-0.2% 

0.3% 

0.1% 

0.3% 

03% 

446.17 

1.2% 

-0.1% 

-0.3% 

0.2% 

0.4% 

0.6% 

453.17 

-0.7% 

-0.1% 

-0.2% 

-0.7% 

-0.6% 

0.5% 


From the above examples it may be seen that the neural network can be used 
as a tool to assess the quality of the individual scans. Particularly, it may indicate 
lack of symmetry in the weld profile, abnormally wide or narrow beads, or other 
undesirable features. Further discussion of neural network techniques may be 
found in [10]. 


The Integrated Weld Quality Evaluation Approach 
The preceding discussion has outlined methods that were developed and 
tested for weld evaluation. For prototype development, the video based system was 
used for the GTAW and GMAW applications. In these, the assessment was based 
upon the quality of the weld beads as determined from the bead images (pixel 
luminosity values and location) alone. The laser scaimer system was used for the 
VPPAW applications. The process evaluation was based upon actual weld 
topography as defined by successive laser scans along the weld axis. 

Specific examples were given in which video images were used to analyze 
welds from the GTAW and GMAW processes and output from a laser scanner was 
used to analyze welds produced by the VPPAW process. As noted in the earlier 
discussion, however, the video-based, the laser-scanning sensing technique, or a 
combination of both may be applied to any of the three processes. 

The laser scanner process can, of course, be applied to the GTAW and 
GMAW weld beads. With appropriate training, a neural network could be applied 
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to the task of feature recognition for these processes as well. By proper treatment 
of the data obtained from the laser scanner, an image of the weld bead may be 
synthesized. This is accomplished by storage, analysis, and display of successive 
laser scans as illustrated in figure 48. This synthetic image may then be analyzed 
and evaluated by numerical algorithms of the type described earlier. 



An overview of the comprehensive weld evaluation system is shown in 
figure 49. In this system, the two components of the overall system are shown 
working in tandem for maximum data acquisition and maximum possible accuracy 
of quality evaluation. However, it has been demonstrated that either side of the 
system will provide a high level of quality of the assessment process. For high 
criticality applications, such as aerospace, both components would be 
recommended. For less demanding applications, such as in consumer products 
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assembly, the video system alone would often be sufficient, reducing overall system 
cost and complexity. 


SYSTEM OVERVIEW 

VIDEO LASER 



Figure 49. An overview of an integrated weld evaluation system. 

This concludes the discussion of weld evaluation approaches applicable for 
the three selected weld processes: GTAW, GMAW and VPPAW. The main 

conclusion from this is that techniques can be developed to assess the quality of the 
weld based on the specified features. Other quality criteria may be built into the 
weld evaluation system if desired, to supplement further the assessment process. 
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CHAPTER VI 

QUALITY CONTROL SYSTEM PROTOTYPE IMPLEMENTATION 

The main components of the prototype system developed for this work are 
outlined in this chapter. Broadly, the system may be separated into hardware 
components, which mostly consisted of commonly available equipment. The 
software, which is the primary deliverable item of this contract, is usable with most 
IBM-PC compatible hardware, as discussed below. This design choice was made to 
encourage commercial replication of the system. 

Hardware Description 

A variety of computer and image processing hardware may be used with the 
software system developed for this work. The primary requirement is that the 
system is based on the IBM-PC architecture and that the computer uses the 80286 
microprocessor or any of its successors. This architecture was chosen to facilitate 
use in commercial applications at modest costs. The specific system components 
which were available at Mid-South Engineering and thus used for this development 
are outlined below. 

The host computer on which the system runs is a Dell System 310 computer, 
which is based on the Intel 80386 miaoprocessor and is equipped with the 80387 
numeric coprocessor. The computer is configured with 4 MB of extended memory 
and a 100 MB hard disk. A Dell VGA computer monitor is used with the computer. 
Its purpose is strictly to control the system through the keyboard or the mouse. 
Graphics display is provided on a separate graphics monitor. The graphics monitor 
is made by Electrohome, but any RGB-compatible monitor would work as well for 
this purpose. 

For image capturing and display the Data Translation DT 2871 Color Frame 
Grabber is used. An optional add-on RGB/NTSC converter, DT 2869, may be used 
to convert images between the RGB and the NTSC television signal standard if an 
NTSC monitor or an NTSC camera is used. The frame grabber receives the red, 
green and blue signals, in addition to the sync signal, from a video camera. 
Internally, the board converts the RGB signal to a hue-intensity-saturation (HSI) 
form, where three 512-times-512 byte buffers store the image information. A fourth 
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auxiliary/overlay buffer is located on the board for data manipulation. The 
DT2871 board outputs the stored image on the RGB format, where the sync signal 
is transmitted with the green-channel or may optionally be obtained through a 
separate connector. 

The video camera used for the work is made by Parasonic. model 
Digital-5000. This camera provides an RGB output signal, in addition to a sync 
.-h^nnel Hs rcsolution is 512-by-512 pixels. The main advantage of this camera 
over many other models, in addition to the high-quality RGB output, is the standard 
35 mm bayonet lens mount. This facilitated the use of eeonomical and widely 
available lenses for this work, including various close-up configurations which were 
necessary for welding inspection. 

A tablet equipped with light sources and a camera column was used as an 
experimental setup to record images of weld samples. Both the light fixtures and 
the camera column were adjustable, permitting the necessary flexibility in adjusting 

lighting conditions and camera position. 

As an alternative or as an addition to inspecting the bead appearance with a 
video camera, a laser profiler may be used, as discussed in context of the VPPAW 
process. Using a laser profiler with the camera gives the system an additional 
degree of reliability, which may be beneficial for exceptionally critical applications. 
A laser-based weld bead profiler is already available for VPPAW work at NASA. 


Software Description 

The software developed under this contract and its operation is described in 
this section. The development of this software constituted a substantial part of the 

work, after the necessary preliminary research. 

The objective of the software was to digitize and process images of welds or 
other inspected material surfaces for the purpose of enhancing discontinuities for 
visual inspection. Furthermore, algorithms were developed to assess the overall 
quality of welds from three distinct welding processes: Gas Metal Arc 

Welding (GMAW), Gas Tungsten Arc Welding (GTAW) and Variable Polarity 
Plasma Arc Welding (VPPAW). The quality is expressed on a numeric scale, e.g., 
from 0 to 10. 

To make the software accessible to a large portion of the commercial user 
market, it was implemented on the IBM-PC platform. Initial work was carried out 
under the MS-DOS operating system. The limitations of early versions of MS-DOS 
were, however, soon encountered. Of particular concern was the limited amount of 
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continuous memory that MS-DOS could access for program data. The MS-DOS 
operating system was limited to a maximum of 640 KB)rtes of memory. This was a 
severe limi tation, as a typical image of 512-by-512 pixels contains a total of 250 
thousand pixels. Multiply this by 3 bytes for each pixel (red, green and blue 
information or, alternatively, hue, intensity and saturation), and the total memory 
requirement for the image is found to be about 750 KB. Not only does the raw 
image exceed the memory space that MS-DOS is designed to handle, but memory 
for the program itself and various data buffers is needed as well. This limitation 
may be circumvented by temporarily directing data to disk or expanded/extended 
memory, but these manipulations call for significant overhead and complications in 
programming. Thus, the OS/2 operating system was selected for the final 
implementation. This operating system is not limited by the older MS-DOS system. 
The downside of choosing a less common operating system, on the other hand, was 
somewhat higher software price and less availability of software tools. Since system 
programming started in the OS/2 environment further developments have taken 
place in the IBM-PC arena. Microsoft Windows, version 3.0, was released, offering 
the menu-driven interface characterizing the OS/2 with its Presentation Manager. 
The MS-DOS system evolved to version 5.0, which greatly enhances memory 
management. As the present work is concluded a new, greatly improved, version of 
OS/2 has been released and Microsoft Windows NT is about to become available. 
Either of these operating system would be considered for further implementation 
and commercial distribution of the image analysis software. 

Menu-Driven Software Overview 

The OS/2 software delivered at the conclusion of this work is described in 
this section. It may be operated on IBM-PC compatible computers, using the Intel 
80286 microprocessor or any of its successors. At least 2.5 MBytes of memory is 
required, but 4 MB is highly recommended. 

The main window of the software system opens up with a choice of five 
general menu items: File, Process, Texture, Edge and Help. By clicking the mouse 
pointer on any of these main selections the user is offered submenus, options or 
other information. The individual items in the menu system will be described in the 
following paragraphs. Note that the monitor displaying video images is referred to 
as the graphics monitor, while the one displaying the menu system and the mouse 
cursor is the system monitor. 
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File: 

Acquire: 

This option allows the user to digitize and store an image from the 
video camera. The computer enables the video camera connected to 
the system and its view is displayed in real-time on the graphics 
monitor. Simultaneously, control buttons are displayed on the system 
monitor which allow the user to cancel the digitization or continue 
and digitize the image on the graphics monitor. The image is assumed 
to be in hue-intensity-saturation form. 

Acquire RGB: 

This option is similar to the one above ( Acquire) , except in this case 
the captured image is in red-green-blue format. This is the default 
operation for image capturing in the present system. 

Open: 

Opens an existing image file for view and further analysis, and 
displays it on the graphics monitor. This is an alternative to the 
acquire commands which import images from the video camera. The 
file that is opened is by default designated by the extension ".IMG". 
The system displays a scroll window from which an existing file (if 
there is any) can be selected and displayed. The file must have been 
acquired earlier and saved by the Save command (see next menu 
item). 

Save As: 

Saves the image that is currently displayed on the graphics monitor. 
This may be an original image, as acquired from the video camera, or 
it may be a processed and altered image. In either case, the system 
asks for a file name to store the image under. Each image file 
requires approximately 787 KBytes of storage space. 

Set Active Buffers: 

This selection prompts the user to select one of the following items: 
Intensity : Selects the memory buffer that contains the image intensity 
information for further processing. 

Saturation : Selects the memory buffer that contains the image 

saturation information for further processing. 
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Hue : Selects the memory buffer that contains the image hue 

information for further processing. 

Overlay : Selects the auxiliary overlay buffer in the image acquisition 
board for further processing. 

Hear Unused : Clears the image buffers that are not currently 

selected for processing. 

Exit: 

Terminates the program. 

Process: 

Blowup: 

This selection prompts the user to select one of the following items: 
Enlarge : Enlarges each pixel on the screen by a factor of two in the 
horizontal and the vertical directions. On first usage, the value of 
every fourth pixel in the image is copied into three neighboring pixels 
adjacent to it, resulting in a formation of squares of four pixels with 
identical values. Successive activations of the enlarging operations 
work in the same manner, creating squares 4-by-4, 8-by-8, etc., of 
identical pixel values. 

Zoom: Similar as the "enlarge" operation, except interpolation is used 
rather than copying of pixel values. 

Integrate: 

This selection prompts the user to select one of the following items: 
Horizontal : Integration of pixel intensities along horizontal pixel 

lines. The result of this operation is a 1-dimensional graph that shows 
the integral of pixel densities as a function of vertical location in the 
image. The user can specify a window of any width (selected as the 
"Number of Lines" option) and horizontal position (selected with a 
scroll bar), extending from the top of the image down to its bottom. 
Each line of pixel within this window is integrated, resulting in a 
measure of the average pixel intensity in that line. This is repeated 
for all lines from the top to the bottom of the image. By clicking the 
mouse on the "Refresh" button in this window the user can display a 
graph of the integrated values. Large values indicate high pixel 
intensities (bright pixels). The displayed graph is normalized so that 
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the largest integral value is plotted at the top of the graph. The 
integration graph may be saved as a disk file by clicking the mouse on 
the "Save" button. 

Vertical : This operation is similar to the one for horizontal 

integration, except integration takes place along vertical lines 
(sometimes referred to as "colunms") in the image, the vertical 
dimension and location of the integration window may be varied, and 
it extends between the left and right boundaries of the image. 

Histogram: 

This selection prompts the user to select one of the following items: 
Intensity : Displays a histogram of pixel intensities in the entire image. 
A low-contrast image would typically have a narrow histogram (not 
much variation in pixel brightnesses) while a high-contrast image 
would result in a broader histogram. The user can select how many 
bins the intensity levels are grouped into. For example, with a total of 
256 intensities offered by the image board, a selection of 64 bms 
would result in pixel intensities 0 through 3 grouped together, 4 
through 7, etc. Clicking the mouse on "Refresh" updates the 
histogram display. 

Saturation : This histogram operation is similar to the one for 

intensity, except a histogram of pixel color saturation levels is 
displayed. 

Hue : This histogram operation is similar to the one for intensity, 
except a histogram of pixel color hue values is displayed. 

Enhance Red: 

Selection of this option converts all hues of red in the graphics 
monitor image into a bright red color. The intensities of all pixels 
with other colors are set to zero, making the effectively black on the 
screen. This operation is useful for analysis of images preprocessed 
with a dye penetrant, which usually is red. 

Equalize: 

The equalization operation generally increases the contrast in the 
image displayed on the graphics monitor. The system accumulates a 
histogram of the original pixel intensity distribution and performs a 
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mapping of it that results in equal distribution of all pixel intensities in 
the image. 

Threshold: 

Pixels with intensities below a specified threshold value will be turned 
off (displayed as black) while those with intensities equal to or above 
this value are displayed as bright white. The threshold value may be 
specified between 0 and 255 by the user through a scroll bar on the 
menu. A click on the "Update" button updates the graphics screen 
image accordingly. 

Average: 

This is a smoothing operation which is useful if the image is "noisy" or 
contains small, irrelevant artifacts. Each pixel intensity on the 
graphics monitor is replaced with the average of the 3-by-3 square of 
pixel of which it is the center. The image on the monitor is updated 
after the averaging. This operation may be iteratively repeated to 
obtain an increasingly smoothed image. 

Inverse: 

The "Inverse" operation replaces each pixel intensity with its 
"opposite" value so that, e.g., all pixels with intensity of 0 become 255, 
1 becomes 254, etc. The effect of this is a "negative" of the original 
picture, similar to what is obtained on negative photographic films. 

Texture: 

Fourier Transform: 

The Fourier transform is a widely used signal processing technique 
which reveals periodic variation in the analyzed signal. The transform 
was implemented in the software as a tool to analyze periodic ripple 
patterns along the weld bead. The result of the Fourier transform of 
the image is a new image, which is displayed on the graphics screen. 
The intensities in the Fourier transform image indicate dominating 
pattern repetitions in the original image. The intensity of the center 
pixel in the Fourier transform is proportional to the average 
brightness of the original image (at pattern fi*equency equal to zero). 
The intensities of Fourier pixels further from the center indicate how 
"dominant" denser pattern repetitions are in the original image. 
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Specifically, a high Fourier intensity at a distance of x pixels to the 
right of the center of the Fourier image, indicates that there is a 
regular horizontal pattern in the original image with a repetition 
period of 512/ (2x) pixels (assuming that the original image has 512 
pixels from left to right). The same holds for analysis of vertical 
patterns. 

8-Rand Fourier Transform : This operation is similar to the "Fourier 
Transform" option, above, except here the image is segmented into 8 
subwindows, each one spanning the full width of the image (from left 
to right) but only l/8th of the total image height. A Fourier 
transform is taken of each of the subwindows. This operation was 
used in early experiments to determine the frequency or spacing of 
weld ripples, but the vertical/horizontal integration procedure was 
later developed and found to be more reliable. 

16 -Rand Fnurip.r Transform : Similar to the "8-Band Fourier 

Transform", except here the image is segmented into 16 horizontal 
subwindows. 

Edge: 

Roberts: 

The edge operations analyze the displayed image and replace all pixel 
intensities with new values that are indicative of brightness variations 
in their neighborhoods. The Roberts operator is essentially the 
gradient operator. Thus, each pixel intensity is replaced by the 
original image gradient at that location. The resulting image 
highlights edges and sharp lines, which characterize many weld 
defects. 

Sobel: 

The Sobel operator is an indicator of the gradients in the image, but 
the gradient is computed from a 3-by-3 neighborhood around each 
pbcel, rather than the 2-by-2 which is used by the Roberts operator. 
The result is an edge finding operator that is less susceptible to 
individual noise pixels than the Roberts operator. 
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Prewitt: 

The Prewitt operator is similar to the Sobel operator in that the local 
gradient is estimated from a 3-by-3 matrix. To improve noise 
performance further, however, it assigns different weights to the pixels 

in the matrix. 

Kirsch: 

Unlike the edge operators listed above, which all estimate the 
gradient by calculating the horizontal derivative and the vertical 
derivative and then combining both, the Kirsch approach also 
determines the two diagonal derivatives before the total gradient is 
calculated. The local matrix used for this calculations is 3-by-3 pixels. 


Help: 

Commands: 

Option for on-line help. 

About Imagine: 

Displays the name of the software and its producer. 

This concludes description of the prototype software system developed for 
this work. In short, in includes many of the tools that provide the foundation for 
visual inspection for quality control. This prototype serves as a template for a 
commercial quality inspection system that is presently under consideration. 
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CHAPTER VII 

CONCLUSIONS AND SUGGESTIONS FOR COMMERCIALIZATION 

The work described in this report has been aimed at exploring means for 
usin g readily available computer technology in the process of quality inspection and 
evaluation. The ffiM-PC computer architecture was selected for the purpose of 
facilitating commercialization for as wide a range of applications as possible. To 
demonstrate the basic principles and potentials of computerized inspection, off-line 
quality control for arc welding was selected as a representative application. Visual 
inspection was emphasized, as it requires relatively modest sensory hardware. 
Furthermore, visual inspection is widely applicable and is in itself an effective 
inspection technique. 

The results of this work illustrate that visual quality inspection may be aided 
and to some extent automated with relatively modest computer resources. A 
prototype software system was developed for the work, which is delivered to NASA 
at the conclusion of this contract. This software may be executed with any 
compatible computer and image processing hardware. Furthermore, a number of 
algorithms and related techniques were developed to aid in evaluating weld quality. 

The next step in this work will be the development of a commercial quality 
inspection system, based on the prototype system developed under this contract. 
The commercial system may be marketed as an integrated package, as well as 
separate hardware and software components. Demonstrations and discussions of 
the prototype system have been conducted with potential commercial customers and 
their coimnents will be used in guiding the development of the commercial system. 
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Historical Overview 

Artificial neural networks are implemented in various configurations and 

some are better suited for certain applications than others. 

One of the earliest accomplishments in neural network implementation was 
the "perceptron" which Rosenblatt designed [11,12]. The perceptron was a primitive 
pattern classifier, designed to recognize optical patterns. A grid of 400 photocells 
(simulating light-sensitive neurons) sensed an illuminated pattern and each 
photocell produced an output signal according to the light intensity it received. An 
array of "associator units" was used to select some of the photocell signals for further 
processing, while others were rejected. The selected signals were transmitted to a 
summing device, and the final output of the perceptron was thus the sum of the 
signals from the associator units. The output of the perceptron was used as an 
indicator of whether the applied optical pattern belonged to the class of patterns 
that the perceptron has been trained to recognize, or not. The associator units, 
which determined which photocell signals were of importance, essentially 
determined the class of patterns recognized by the perceptron. The perceptron had 
certain limitations, the most severe one being that it could not learn to classify any 
arbitrary patterns. This lack of generality was overcome later by the introduction of 
multiple-network structures and adding a nonlinearity transfer function to the 
summing node. 

Minsky and Papert analyzed the basic perceptron further [13] and 
demonstrated that, for example, that perceptrons could not distinguish connected 
pattern figures from unconnected ones. Furthermore, they showed that the basic 
perceptron can not compute functions that are not linearly separable, such as the 
Exclusive-OR Boolean function. Multiple-layer perceptrons, using nonlinear 
transfer functions, were introduced to overcome the limitations of the basic 
perceptron. Although these structures were shown to be more powerful than the 
original perceptron, they did not gain general interest until the backpropagation 
algorithm was devised [14,15,16]. This algorithm provided means for systematically 
adjusting the various weights of the multi-layer network until it had been adequately 
trained. 
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Figure 50. The processing element of a neural network. 

The multi-layer networks consist of a number of processing elements, which to 
a certain degree resemble the connected neurons in biological systems (refer to 
Figure 50). A processing element accepts one or more signals, which may be 
produced by other processing elements or applied externally (e.g., provided by a 
process sensor). The various signals are individually amplified, or weighted, and 
then sununed together within the processing element. The resulting sum is applied 
to a specified transfer function, and the function value becomes the output of the 
processing element. 


The Backpropagation Network 

A neural network and its adaptation procedure using the backpropagation 
algorithm is best illustrated by an example. 
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Figure 51. A backpropagation neural network. 

Figure 51 shows a small neural network consisting of 8 nodes, arranged in 
two hidden layers of 3 nodes each and one output layer of 2 nodes. Each node i in 
the first hidden layer produces a single numeric output which we denote as xf^K 
Similarly the nodes of the second hidden layer are labelled through x^^K The 
3 inputs and 2 outputs of the network are xq through X 2 and yQ through yj 
respectively. Each node accepts numeric data through a number of input links, each 
of which multiplies the input data with a weight factor. The weight factor associated 
with the link from to the node producing is aimotated as and a similar 
convention holds for the links between other layers. Each node calculates its output 
by summing its weighted inputs and using the result s as the argument of a nonlinear 
function associated with the node. For our application this function is the same for 
all nodes: 

f(s) = [1 + exp[-(s-c)]]-^ 


( 11 ) 
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where s is the sum of the node inputs and c is an internal offset value. Qearly the 
node output will be confined to the range 0 < f(s) < 1. Because the limiting values, 
0 and 1, will only be approached as s approaches +/- infinity all input and output 
data are scaled so that they are confined to a subinterval of [0...1]. A practical 
region for the data is chosen to be [0.1. ..0.9]. In this case each input or output 
parameter p is normalized as Pn before being applied to the neural network 
according to: 

p. = [(0.9-0.1)/(p„^p^)](p-p„i„)+0.1 (12) 

where Pn,ax ^^d p^in are the maximum and minimum values, respectively, of data 
parameter p. The network starts calculating its output values by passing the 
weighted inputs to the nodes in the first layer. The resulting node outputs of that 
layer are passed on, through a new set of weights, to the second layer, and so on 
until the nodes of the output layer compute the final outputs. 

Before practical application, the network has to be trained to perform the 
mapping of the three input parameters to the two output parameters. This is done 
by repeatedly applying training data to its inputs, calculating the corresponding 
outputs by the network, comparing them to the desired outputs, and altering the 
internal parameters of the network for the next round. The training starts by 
assigning small random values to all weights (wjj) and node offsets (c^) in the 
network. The first three input data values are presented to the network which in 
turn calculates the two output values. Because the initial weights and node offsets 
are random these values will generally be quite different from the desired output 
values, Dq and Dj. Therefore the differences between the desired and calculated 
outputs have to be utilized to dictate improved network values, tuning each weight 
and offset parameter through back propagation. The weights preceding each output 
node are updated according to 

Wij(t+1) = wjj(t) + ndjx/2) (13) 


where n is a correction gain and dj is the conection factor 
dj = yj (1-yj) (dj-yj) 


( 14 ) 
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Oearly, each weight will be increased if the calculated output from its node 
is less than the desired value, and vice versa. The correction factors used to update 
weights preceding hidden layer nodes are updated according to 

dj = Xj (l-Xj)2 (dk-Wjk) (15) 

k 

where the k applies to the node layer succeeding the one currently being updated. 
The offset parameter c of each node is treated as an additional weight factor and 
updated in the same manner. 

The weights and offsets of the neural network are recalculated during the 
backpropagation as outlined above. Then the network repeats calculation of output 
values based on the same input data, compares them to the desired output values, 
and readjusts the network parameters through yet another backpropagation phase. 
This cycle is repeated until the calculated outputs have converged sufficiently close 
to the desired outputs or an iteration limit has been reached. Once the neural 
network has been tuned to the first set of input/ output data, additional data sets can 
be used for further training in the same way. To ensure concurrent network 
adaptation to all sets of data the entire training process may be repeated until all 
data transformations are adequately modeled by the network. This requires, of 
course, that all the data sets were obtained from the same physical process and 
therefore the underlying input/output transformation is consistent. 

As noted above, the training iteration process may be terminated either by a 
convergence limit or simply by limiting the total number of iterations. In the former 
case we use an error measure e defined as following; 


e 




M-l 

{ ^ (dk,m"yk,m) } 

m=0 


( 16 ) 


where K is the number of input/output data sets used for training, M is the number 
of network output parameters in each data set, and (d^ niyk,m) is the error in the 
network calculation of parameter m in data set k. The error measure, e, changes 
after each round of network weight adjustments. In the long run e decreases as the 
network is refined by training iterations. Using this indicator one can program the 
network to terminate the iterative tuning process as soon as e reaches some 
threshold value, eQ. Alternatively, a given network may not be able to reduce the 
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error measure down to the specified eo- In that case the iterations may be 
terminated by simply specifying a maximum number for them. 

The training mode, as described above, is a precondition for actually applying 
the neural network in the application mode. In this mode entirely new input data is 
presented to the network which, in turn, predicts new outputs based on the transfer 
characteristics learned during the training. If this new data is obtained from the 
same local region of operation of the process as during the training phase, data from 
the input/output relations should be governed by the same underl)dng process and 
the neural network should perform adequately. The neural network is not updated 
in the application mode. 

The above discussion should suffice to give an overview of neural networks 
and how they work. The reader is referred to the bibliography in this report for 
further information on the mathematics and theory of artificial neural networks. 

Neural Network Applications 

Artificial neural networks have been used in a number of applications, 
ranging from radar signal processing to financial analysis to robotics. 

One of the earliest application of neural networks is pattern classification. 
Several military applications are currently under development, where neural 
networks are implemented to recognize targets in real-time. It is believed that a 
proximity fuze for a smart weapon can be realized by a multi-layer neural network, 
which in turn, can be implemented on a single silicon chip [17]. 

Process modeling and control, is one are in which neural networks have 
proven applicable. The application of neural networks for modeling and control of 
the arc welding processes has been pioneered by the authors of this report and the 
results have been published in [18,19,20] 

Neural networks are likely to appear in robotic manipulators for trajectory 
control. Simulations have proven to be successful and it appears to be only a matter 
time before specific hardware becomes available for this purpose. 

Finally, signal processing is one are where neural networks have been tried 
with some success. Because neural networks provide an adaptive, nonlinear 
capability, they appear to offer a promising alternative to more traditional efforts in 
this area. 
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Neural Network Hardware and Speed Considerations 

A number of hardware accelerators and dedicated neural network computers 
are presently in the design phase or prototype production stage. Furthermore, a 
number of neural network software vendors have optimized their code for generic 
numerical accelerator boards or chips. The use of such dedicated hardware 
components may be desirable if substantial amount of re-training is required in a 

production system. 

NeuralWare, Inc. [21] developed the software that was used for most of the 
neural network simulations presented in this work. This vendor offers a genenc 
numerical accelerator board (ALACRON) which increases the execution rate of the 
neural networks by a factor of approximately an order of magnitude. Using this 
board, typical execution times for a 15-node network with 2 inputs and 4 outputs is 
less than 1 millisecond. 

More experimental neural network hardware includes the Network 
Emulation Processor, developed by Claude Cruz-Young, who was previously at IBM 
Research Laboratories. This is a medium-scale VLSI-based virtual processor. It is 
currently used as an emulator for the Parallel Associative Network (PAN) 
workstation with an IBM PC. 

ANZA and ANZA Plus are two boards designed specifically for neural 
network simulation, and developed at the Hecht-Nielsen Neuro-Computing 
Corporation. Both boards are based on the Motorola 68xxx series processors and 
co-processors. The ANZA board can achieve up to 45,000 inter-connects per 
second, while the ANZA Plus reaches as high as 1,500,000 inter-connects per 
second. 

Neural networks are, although on a small scale so far, implemented on 
dedicated silicon chips. AT&T has been developing a neural computing chip that 
has 256 neuron-transistors and more than 100,000 synapse-resistors. These resistors 
are "burned in" at the time of manufacturing, so training and simulation of the 
application must be completed before the chips are ordered. AT&T is currently 
working on a re-programmable chip. 

It is important to distinguish between the training speed and execution speed 
of neural networks. The training speed of neural networks is generally very low. 
During the training process the internal parameters of the network are iteratively 
adjusted until the network training has been optimized. For the backpropagation 
network the network weights are the adjustable parameters. The training is usually 
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carried out off-line and optimized before the network is finally applied in the 
application. The task of weld profile analysis with neural networks, which was 
demonstrated for this work, was found to be executed in a fi*action of a second on 
relatively modest computer systems, which is adequate for most inspection 
applications. 



APPENDIX B 


Source Code Listing for the 
Quality Analysis System 
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# File HKJMAGE 

CFLAGS = -c -AL -G2sw -Owx -W3 -FPi87 

CFLAGSC = -c -AL -G2s -Ox -W3 -FPi87 

#CFLA6S = -c -AL -G2sw -Zi -Od -W3 -FPi87 

#CFLAGSC » -c -AL -G2s -Zi -Od -U3 -FPi87 

all: iinage.exe 

image. obj: image. c dialog. h menu.h 
cl $(CFLAGS) image. c 

process. obj: process. c 

cl $(CFLAGSC) process. c 

edge. obj; edge.c 

cl S(CFLAGSC) edge.c 

integ.obj: integ.c 

cl S(CFLAGSC) integ.c 

hi St. obj: hist.c dialog. h 
cl $(CFLAGS) hist.c 

fourn.obj; fourn.c 

cl S(CFLAGSC) fourn.c 

texture. obj: texture. c 

cl $<CFLAGSC) texture. c 

image. res: image. rc dialog. h image. ico menu.h dialog. dig 
rc -r image 

image.exe: image. obj process. obj edge. obj hist. obj mapn.obj copy. obj image. def integ.obj fourn.obj 

texture. obj 

link /CO /align: 16 Simage.rsp 
rc image. res 

image.exe: image. res 
rc image. res 
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/* 

File DIALOG. H • 


#6ef ine 

IDB^AUTO 

1013 

#def ine 

lOB^SAVE 

1012 

#def ine 

lOT^INT 

1011 

#def inc 

IDEJNT 

1010 

#def ine 

IDSJNT 

1009 

#def ine 

IDB.CANCEL 

2 

#def ine 

1DB_0IC 

1 

#def1ne 

ID^RESOURCE 

999 

#def ine 

ID^HISTO 

120 

#def ine 

lOO^ABCXJT 

1001 

#def ine 

IDO^SCROLL 

1005 

#define 

IDD.BUF 

1002 

#def ine 

lOO^HISTO 

1003 

#def ine 

IDD.INTEGRATE 

1004 

#def ine 

IDD_PATH 

110 

#def ine 

IDD^FILEEDIT 

111 

#define 

IDD^DIRLIST 

112 

#define 

IDD_FILELIST 

113 

#define 

IDD^HEAD 

114 

#define 

IDD^SCROLLHEAO 

130 

#def ine 

IDD^SCROLLBAR 

131 

#def i ne 

IDD.SCROUNUH 

132 

#def ine 

I DC J NT 

140 

#def ine 

IDC^SAT 

141 

#def ine 

IDC^HUE 

142 

#define 

IDC^OVL 

143 

#def ine 

IDC^CLR 

144 

#def ine 

IDW_HUINI 

150 

#def ine 

IDU.HUINH 

151 

#def ine 

IDU^HUINS 

152 

#def ine 

IDB_16 

153 

#def ine 

I0B_32 

154 

#def ine 

IDB_64 

155 

#def ine 

I OBJ 28 

156 

#def i ne 

IDB_256 

157 

#def ine 

IDBJUE 

158 

#def ine 

IDBJAT 

159 

#def ine 

IDB^INT 

160 

#def ine 

IDDJILEIO 

1000 

#def ine 

IDUJAR 

1006 

#def ine 

IDB_UPDATE 

1008 
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/• File EDGE.H */ 

int sobel_inask[8] [31 [3] -{ 

{C2.1,0).{1,0,-1>,{0,-1,-2». 

«1.0,-1>,{2,0,-2>.<:i,0,-1», 

« 0 ,- 1 ,- 2 },{ 1 , 0 .- 1 >,{ 2 , 1 , 0 », 

«- 1 ,- 2 .- 1 >.< 0 , 0 , 0 >,{ 1 . 2 . 1 ». 

«- 2 ,- 1 . 0 },<- 1 , 0 , 1 >,{ 0 , 1 , 2 », 

«-1,0,n,{-2,0.2>,{-1,0.1», 

« 0 , 1 , 2 }.{- 1 , 0 , 1 >,{- 2 ,- 1 , 0 » 

>; 

int prewitt_Mask[8] [3] [3] = { 

>; 

int kirschjnask[81 [31 [3] s { 

{{5.5,5>,{-3.0,-3>,{-3,-3.-3», 
«5,5,-3>,<5,0,-3>,{-3,-3,-3)>. 
«5.-3,-3>.<5.0.-3>.{5,-3.-3)>, 
«-3.-3,-3>,{5,0.-3),{5,5,-3», 
«-3.-3,-3>,{-3.0,-3),{5.5,5)). 
{{-3, -3, -3>, {-3, 0,5), (-3, 5,5)}, 
{(-3, -3,5), {-3, 0,5), (-3, -3,5)), 
{{-3,5, 5), {-3, 0,5), {-3, -3,-3)) 
); 



86 


/* File MENU.H V 

#define IDH_FILE 10 
#define IDM_ACQUIRE 11 
«define IDM_ACQUIRERGB 12 
#define IDN_OPEN 13 
«define IDN_SAVEAS U 
#define IDN_BUFFERS 15 
#define I0M_EXIT 16 

«def1ne IDM_PROCESS 20 
«define IDM_BLOUUP 21 
«define IOM_ZOOM 22 
«define IDM_ENLARGE 23 
#define IDH_ENHRED 24 
Bdefine IDM_EQUALIZE 25 
Bdefine IDM_THRESHOLD 26 
Bdefine IDM_AVERAGE 27 
Bdefine IDHJNTEGRATE 28 
Bdefine IDH_HOR 29 
Bdefine IDM_VER 30 
Bdefine IDM_HISTO 31 
Bdefine IDM_INT 32 
Bdefine IDM_SAT 33 
Bdefine IDM_HUE 34 
Bdefine IDM_IHVERSE 35 

Bdefine IDM_EDGE 40 
Bdefine IDH_ROBERTS 41 
Bdefine IDN_SOBEL 42 
Bdefine IDH_PREWITT 43 
Bdefine IDH_KIRSCH 44 

Bdefine IDH_TEXTURE 50 
Bdefine IDM_ONEITER 51 
Bdefine IDM_HULT I THRESH 52 
Bdefine IDH_DOITALL 53 
Bdefine IDM_FOURIER 54 
Bdefine IDH_BAND_8 55 
Bdefine IDH_BAND_16 56 

Bdefine IDH_HELP 60 
Bdefine IDN_HELP_COMMANDS 61 
Bdefine IDN ABOUT 62 


Bdefine ID FRAME 99 
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/* File COPY.C V 

void copy(unsigned char huge unsigned char huge *); 

void copy(ipic,opic) 
unsigned char huge *ipic; 
unsigned char huge *opic; 

unsigned int i; 
for (isO;i<4;i*M-) i 
opic 65536L; 
ipic +s 65536L; 

*opic « *ipic; 

> 

> 
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/• File EDGE.C */ 

#include <math.h> 

#include “util.h” 

#includc “edgc-h** 

#define HEIGHT 512 
#define WIDTH 512 

extern unsigned char huge ipiclHEIGHT) EUIDTH3 ; 
extern unsigned char huge opicCHEIGHT] [WIDTH]; 

extern int ha; 
void robert8_main(void); 
void sobel_iRain(void); 
void prewitt_main(void); 
void kirsch_main(void); 

unsigned char mask_max(int <*) C3] [3] ,int, int); 
void maskit(int mask [8] [3] [3]); 

void robert8_main(void) < 

int i,j,d1,d2,p; 

unsigned char huge *tpic; 

tpic = (unsigned char huge *) ipic; 

au_readscr(ha, ipic,0); 

for Ci*0; i<HEIGHT-1; i++.tpic-»'+) ( 

for (j=0; j<UI0TH-1; j++,tpic-M-) { 

d1 » (int) *(tpic+1) - (int) *tpic; 
d2 = (int) *(tpic+WIDTH) - (int) *tpic; 
p = (int) (.707 * sqrt((double) d1*d1+d2*d2)); 
♦tpic * (unsigned char) p; 

> 

> 

au_wri tescr(ha, ipic,0); 

> 

void maskit(mask) 
int mask [8] [3] [3] ; 
i 

int i,j; 
int a,c,t; 
int (*m)(3); 

unsigned char huge *topic = (unsigned char huge *) (opic[1]); 
au_readscr(ha, ipic,0); 
for (i=1; i<HEIGHT-1; i«M., topic+=2) ( 
for (j=1; ]<W10TH-1; topi C++) ( 
for (t=0,a=0;a<4;a++) ( 
in s mask [a] ; 

c * ((int)ipic[i-1] [j-1] )*m(0) [0] ; 
c += ((int)ipic[i] (j-13 )*mC0] [11 ; 
c +» ((int)ipic[i+1) [j-13 )*m[0] [23 ; 
c +* ((int)ipic[i-13 [jl )*m[13 [03 ; 
c += ((int)ipic[i3 [jl )*m[13 [13 ; 
c +* ((int)ipic[i+13 [j3 )*m[l3 [23 ; 


^'- 3 . 


c ((int)ipicCi*1] [j+1] >*mE2] [01 ; 
c ((int)ipictil [j>1] )*m[2] Cl] ; 
c ♦= «int)ipic[i*11 [j+1])*m[2] [21; 
c * abs(c)/A; 
if (ot) t*c; 

> 

•topic ■ (unsigned char) t; 

> 

> 

au_writescr(ha,opic,0); 

> 

void 8obel_fnain(void) i 
maski t ( sobe l_inaslc) ; 

> 

void prewi tt_main(void) { 
mask 1 t ( preu i t t_mask ) ; 

> 

void kirsch_main(void) < 
mask i t ( k i rsch_mask ) ; 

> 
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/* File FOURN.C V 

^define INCL^WIN 
#define INCL.DOS 
#include <os2.h> 

#include <math.h> 

#include <util.h> 

#define SWAP(a,b) tempr*(a);(a)=(b);(b>*tenpr 
void fournCfloat huge *,long *,int,int); 
extern HWNO hwndFrame; 
extern unsigned char huge tpic[512] [512] ; 
extern unsigned char huge opic[512] [512] ; 
extern int ha; 

void fft(void) 
i 

float huge *arr; 

float huge *tarr; 

float tf 1,tf2,mx=(f loat)0.; 

SEL selHuge; 

static long dims[2]K512L,512L>; 
long li,lj; 

unsigned char huge *tp = (unsigned char huge *)tpic; 

if (DosAllocHuge(32.0,&selHuge,0,SEG_N0NSHARE0)) i ; /* allocate 2Meg of memory V 
DosBeepC 2000 , 1 000 ) ; 

Wir#lessageBox(HWND_0ESKT0P, hwndFrame, "Could not allocate memory", NULL,0,MB_OK | 
HB^ICONEXCLAMATION); 
return; 

> 

tarr ■ arr * (float huge ♦) MAKEP(selHugc,0); 

0osBeep( 8000 , 300 ) ; 

au_readscr(ha,tpic,0); /* read intensity data V 
for (li*0L;li<262144L;li++) < 

♦tarr * (float) *tp; 

tarr-M-; 

♦tarr = (float)O.; 

tp^^; 

tarr+-»*; 

> 

tarr = arr; 

tp = (unsigned char huge ♦)tpic; 
fourn(arr,dims,2,1); 
for (li=0L;li<262144L;li-**+,tarr+=2) i 
tfl = ♦tarr; 
tf2 * ♦(tarr+1); 

♦tarr * (float) log(sqrt(tf1^tf 1+tf2^tf2)); 
if (♦tarr > mx) mx = ♦tarr; 

> 

tarr * arr; 

for (li*0L;li<262144L;li+"*‘,tarr^=2,tp<^) { 
tfl * ♦tarr/mx^(f loat)255.; 
if (tfl < (float)O.) tfl * (float)O.; 
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*tp * (urtsigned char) tf1; 

> 

for <li»0L;li<256L;li^) 

for (I j«0L;l j<256L;l < 

opictli+256L]Clj+256L] = tpictli] tl j] ; 
opictliHUJ * tpicCli+256LHlj+256L]; 
opictli+256Ll[ljJ * tpicm]Clj+256L]; 
op?c[lil tlj>256L3 « tpicCli+256L) [I j]; 

> 

au>ritcscr(ha,op1c,0); write intensity data */ 

DosFreeSeg(selHuge); 

> 

void bfft(bands) 
int bands; 

< 

float huge *ar re- 
float huge *tarr; 
float tf1,tf2,inx*(float)0.; 

SEL selHuge; 

static long dimsI23*<64L,512L>; 
long li,lj,lk; 

unsigned char huge *tp * (unsigned char huge *)tpic; 
dims CO] » 512L/(long)bands; 

if (Do8AllocHuge(32.0,»selHuge,0,SEG_NONSHARE0)) < ; /* allocate 2Heg of men»ry V 

DosBeepC 2000 , 1 000 ) ; 

Uirtle88ageBox(HUND_DESICT0P,hwndFrame, "Could not allocate aie(iiory",NULL,0,MB_OIC | 
MB_I CONEXCLAMAT I ON ) ; 
return; 

> 

tarr » arr ■ (float huge *> MAKEP(selHuge.O); 

DosBeep( 8000 , 300 ) ; 

au_readscr(ha,tpic,0); /* read intensity data V 
for ( I i«0L; I i<262144L; li+'*‘> i 
*tarr * (float) *tp; 
tarr++; 

♦tarr * (float)O.; 

tp*+; 

tarrM-; 

> 

tp * (unsigned char huge *)tpic; 
for {li*0L;li<512/dimsC03;li‘*^) i 

tarr * &arr [li*dimsC03*dimsCl3*2] ; 
fourn(tarr,dims,2,l); 

> 

tarr = arr; 

for (li*0L;li<262144L;li++,tarr+=2) i 
tf1 = *tarr; 
tf2 = *(tarr+1); 

♦tarr « (float) sqrt(tf1*tf Utf2^tf2); 
if (♦tarr > (float)l-) 

♦tarr * (f loat)log(*tarr); 

else 


♦tarr * (float)O,; 
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if (*tarr > nw) mx * *tarr; 

> 

tarr ■ arr; 

for ( li*0L; li<26214AL; li++,tarr+=2,tp^+) { 
tf1 * *tarr/mx*(f loat)255.; 

*tp = (unsigned char) tfl; 

> 

for (lk*0L;lk<512/dims[0];lk++) i 
for (li=0L;li<dimsC01/2;li++) i 
for (lj*0L;lj<256L;lj++) i 

opicUk^dimsCOl+li] Cl j] * tpicClk*dims[03+diinsC0J/2+li] tlj+2563; 
opicClk*dims[0]^dimst0]/2+li3 Clj+2563 * tpicClk*dims[03 + l i3 Clj3; 
opic[lk*dimsC03 + li3 Clj+256] = tpicClk*dinsC03+dimsC03/2+li3 Clj3; 
opic[lk*dimsC03+dimst03/2+li3 [lj3 * tpic[lk*dimsC03 + U3 Cl j+2563 ; 
> 

> 

> 

au_writescr(ha,opic,0>; /* write intensity data V 
DosFreeSeg(selHuge); 

> 


void fourn(data,nn,ndim, isign) 

float huge *data; 

int ndim, isign; 

long *nn; 

i 

long i 1 , i 2 , i 3 , i 2rcv, i3rcv, i p1 , i p2 , i p3 , i f pi , i f p2; 
long ibit,k1,k2,n,nprev,nrem,ntot; 
int idim; 

float tenpifternpr; 

double theta, wi,wpi,wpr,wr, wtefnpf- 

ntotsiL; 

for (idim=0;idim<ndim;idirrH-+) 
ntot *= nnCidim]; 
nprev=1L; 

for (idim=ndim-1; idim>=0; idim--) { 

DosBeep( 1 000 , 500 ) ; 
n=nn[idim] ; 
nrenpnt ot/ ( n*nprev ) ; 
ip1*nprev « 1; 
ip2=ip1*n; 
ip3=ip2*nrem; 
i2rev=1; 

for (i2=1; i2<=ip2;i2^=ip1) { 
if (i2 < i2rev) ( 

for <i1=i2; i1<=i2+ip1-2;i1>=2) < 

for (i3=i1;i3<=ip3;i3>=ip2) i 
i3rev=i2rev+i3- i2; 
SUAP(data[i3-l3 ,dataCi3rev-1] ); 
SWAP(data[i33 ,dataCi3rev] ); 

> 


> 
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> 

ibit»ip2 » 1; 

while (ibit >* ip1 && i2rev > ibit) { 
i2rev -= ibit; 
ibit »- 1; 

> 

i2rev +« ibit; 

> 

ifplaipl; 

while <ifp1 < ip2) < 

ifp2=ifp1 « 1; 

theta= i s i gn*6 . 2831853071 7959/ ( i f p2/ i pi > ; 

wt enrp=s i n( 0.5* theta); 

wpr * -2.0*wtefnp*wtenp; 

wpiBsin(theta); 

wr*1 .0; 

wisO.O; 

for (i3=1;i3<=ifp1;i3*=ip1) < 

for (i1*i3;i1<=i3+ip1-2;iU=2) ( 

for (i2=i1;i2<=ip3;i2+=ifp2) i 
k1=i2; 
k2=kUifp1; 

tenpr=wr*data[k2-11-wi*data[k23 ; 
tenpi swr*data [k2] +wi *dat a [k2- 1 3 ; 
data [k2- 13 =data Ck1 - 13 - tempr; 
data Ik23 =da ta Ck1 3 - tenpi ; 
datatk1-13 *= tenpr; 
dataCk13 tenpi; 

> 

> 

wr=(wtenp=wr)*wpr-wi*wpi+wr; 

wi=wi*wpr+wtenp*wpi+wi; 

> 

ifp1=ifp2; 

> 

fprev *= n; 

> 

> 


#undef SWAP 
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/♦ File HIST.C V 

#define INCL^UIM 
#define INCL^DOSPROCESS 
^include <os2.h> 

#inciude <stdio.h> 

#include **dialog.h” 

extern HUND hwndFrame; 
static char buf[80]; 

MRESULT EXPENTRY HistoWndProc(HWNO hwnd, USHORT msg,MPARAM mp1,HPARAM mp2) 

< 

HP5 hps; 

USHORT i; 

ULOMG width, height; 
float incxjfx; 

POINTL pt; 

SWP swp; 

unsigned int num_hist_ levels; 
float *hist; 
switch(msg) i 

case UM.PAINT: 

mjit_h is thieve Is * WinOueryWindo«UShort(hwnd,8); 
hist > (float *) UinQueryUindowPtr(hwnd,0); 
hps s WinBeginPaintChwnd, NULL, NULL); 

GpiErase(hps); 

GpiSetPattern(hps,PATSYM_0IAG1); 

Wi nOueryWi ndowPos ( hwnd, &sup) ; 
width s (ULONG) swp.cx; 
height » (ULONG) swp.cy; 

incx z (float) width / (float) num_hist_levels; 
pt.x s OL; pt.y=0L; 

Gpi Hove( hps , &pt ) ; 

pt.x * width-IL; pt.y = height-IL; 
GpiBox(hps,DRO.OUTLlNE,&pt,Ol,OL); 
fx = (float) 0.; 
pt.x = OL; 

for (isO;i<nLin_hist_levels;i'»-*‘) { 
pt.y = OL; 

Gpi Hove ( hps , &pt ) ; 
fx ♦= incx; 
pt.x = (long) fx; 

pt.y ♦= (long) ( (float) height ♦ histCi]); 

Gpi Box (hps , DROPOUT L I NE , Apt , OL , OL ) ; 

DosBeep( (USHORT) ((f loat)height*(f loat)10.*hist li) ), 10); 

> 

WinEndPaint(hps); 
return 0; 

case UM_ERASEBACXGROUND: 
return 1; 

> 

return UinOefUindowProc(hwnd,msg,mp1,inp2); 
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> 


MRESULT EXPENTRY IntWndProcCHWND hwnd, USHORT msg,MPARAM mp1,MPARAM mp2) 
i 

HPS hps; 

USHORT i; 

ULONG width, height; 
float incXffx; 

POINTL pt; 

SUP swp; 
float *inte; 
switch(msg) < 

case UH^PAINT: 

inte * (float •) UinQueryWiridowPtr(hwrxJ,0); 
hps > UinBeginPaintfhund, NULL, NULL); 

GpiErase(hps); 

Gpi SetPat tern( hps , PATSYM^D lAGl); 

U i nQueryW i ndowPos ( h wnd , &s wp) ; 

width s (ULONG) swp.cx; 

height ^ (ULONG) swp.cy; 

incx = (float) width / (float) 256; 

pt.x = OL; pt.y=0L; 

GpiMove(hps,&pt); 

pt.x = width-IL; pt.y = height-IL; 

Gpi Box(hps ,DRO_OUTL I NE ,ipt , OL , OL ) ; 
fx = (float) 0.; 
pt.x = OL; 

for (i=0;i<256;i'*'+) < 
pt.y * OL; 

Gpi Move( hps , Apt ) ; 
fx incx; 
pt.x = (long) fx; 

pt.y ♦= (long) ( (float) height * inte(i)); 

Gpi Box( hps , DROPOUT L ] NE , &pt , OL , OL ) ; 

> 

UinEndPaint(hps); 
return 0; 

case UM.ERASEBACKGROUND: 
return 1; 

> 

return WinDefWindowProc(hwnd,msg,nip1 ,irp2); 

> 
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/* File IMAGE. C V 

«define INCL.WIN 
#include <os2.h> 

#include <stdio.h> 

#includc <string.h> 

#include **dialog.h" 

#include “roenu.h” 

#include "procdef.h" 

#include "util.h** 

void far _loadds PMcrr(char *); 

void enlarge_main(void); 

void zoom_iBain(void); 

void enhred_mBin(void); 

void equal ize_niain( void); 

void roberts_main(void); 

void 8obel_main(void); 

void prewi tt_main<void); 

void kirsch_main(void); 

void init_thresh(unsigr>ed char); 

void cal I ^thresh (unsigned char); 

void end_thresh(int); 

void aver age( void); 

void inverse(void); 

void texaverage(void); 

void calc_hist(f loat unsigned int, unsigned int); 
void integ(USKORT,USHOftT,USHORT); 
void fft(void); 
void bfft(int); 


extern float histi (2561 ,hists [256] ,histh [2561 ; 
extern float i nteg_h (2561 , int eg^v [2561 ; 

unsigned int act_buf=0x07; /* bits indicate active buffers V 

#def ine BUSY WinSetPointer(HWMD_DESICTOP,WinQuerySysPointer<HWNO_DESKTOP,SPTR_UAIT, FALSE)) 

#def i ne NOTBUSY WinSetPointer C HWND_OESKTOP, Wi nQuerySysPointcr(HWND_DESKTOP , SPTR_ARROW, FALSE ) ); 

#define BNOTB(x) BUSY; x; NOTBUSY 

#define max(a,b) ( (a) < <b) ) ? (b) : (a) 

#define min(a,b) ( (a) < (b) ) ? (a) : (b) 

SHORT ParseFileName (CHAR CHAR *) ; 

char szFileName [801 ; 

int ha; 

int thresh=128; 
char buf[51; 

MRESULT EXPEMTRY ClientWndProc(HWND,USHORT,MPARAM,MPARAM); 

MRESULT EXPEMTRY HistoWndProc(HWND,USHORT,MPARAM,MPARAM); 

MRESULT EXPEMTRY lntWndProc(HWND,USHORT,MPARAM,MPARAM); 


MRESULT EXPENTRY AboutDlgProc(HUND,USHORT,HPARAM,HPARAM); 

MRESULT EXPENTRY Operi)lgProc(KUND,USHORT,MPARAH,MPARAM); 

MRESULT EXPENTRY ThreshDlgProc(HWND,USHORT,MPARAM,MPARAM); 

MRESULT EXPENTRY BufDlgProc(HWND,USHORT,MPARAM,MPARAM); 

MRESULT EXPENTRY lntDlgProc(HUND,USHORT,MPARAM,MPARAM); 

MRESULT EXPENTRY HstDlgProc(HWND,USHORT,MPARAM,MPARAM); 

HWND hwndFrame,hwndClient; 

HAB hab; 

HMQ hmq; 

QMSG (|nsg; 

static CHAR szClientClassns"Iniage*'; 

static CHAR szHistoClassC]s*'HistoClass'*; 

static CHAR szIntClassHs^lntClass*'; 

static CHAR title[]=**IMAGine -- Inspection Software”; 

static USHORT tefnp; 

int main(void) { 

hab=UinInitiaUze(0}; 
hmq=Wi nCreateMsgOueueC hab, 0 } ; 

U i nReg i sterC I ass ( hab, sz I nt C I ass , I ntUndP roc , OL , 
sizeofCfloat *)+sizeof<USHORT)); 

Ui nReg i sterC I ass< hab, szC I i entC I ass , C I i entUndProc , OL , 0) ; 

Ui nReg i sterC I ass( hab, szH i stoC I ass , H i s toUndProc , OL , 
sizeof (float *)^3*sizeof (USHORT)); 

hwndFrame ^ Wi nLoadO I g( KUND_.DESI(TOP , HUND^DESXTOP , NULL , ( HMOOULE )NULL , I D_FRAME , NUL L ) 

UinLoacMenu(hwndFrame,(HM(X>ULE)NULL,ID_RESOURCE); 

UinPostMsg(hwndFraine,WH_UPDATEFRAHE,HPFR0H2SH0RT(FCF_MENU,0),NULL); 

Ui nSetUi ndowT ext (hwndF rame, title); 

UinSetFocus(HWND_DESXTOP,WinUindouFronilD(hwndFrame,FID_CLIENT)); 

au_init(&ha); 

au_crr_msgs(ha,AU_ERR_QUIT|AU_ERR_PRINT,PMerr); 
au_out put ( ha , 0x8000 ) ; 

whi le (W{nGetMsg(hab,&qmsg, NULL, 0,0)) 

Win0ispatchMsg(hab,&<7nsg); 

U i nOes t royU i ndou( h undF rame ) ; 

U i nO es t royMsgQueue ( hmq ) ; 

UinTerminate(hab); 
au_out put ( ha , 0x0000 ) ; 
au_end(ha); 
return 0; 

> 


MRESULT EXPENTRY ClientUndProc(HWND hwnd, USHORT msg,MPARAH mp1,MPARAH mp2) 
suitch(msg) i 
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case WH^CREATE: 

Wi nPostMsg(Ui nUi ndowF romlD ( Kund, IDB_64 ) , BH_SETCHECK , 

HP FR0M2SH0RT ( TRUE , 0 ) , NULL ) ; 
return 0; 
case UM_C0HHAND: 

8witch(C0MHANDMSG(&msg)->and) i 
case IDM^EXIT: 

U1 nPos tMsg( hwnd , UM^QU I T , OL , OL ) ; 
return Gr- 
ease IDM^HUE: 
case IDM^SAT: 
case IDMJNT: 

temp - COMMANDMSG(&msg)->cnd - IDH^INT; 

WinShowUinck>w(WinLoadDlg<HWND_OES»CTOP,hwnd,HstDlgProc,(HMOOULE)NULL. 

lOO^HISTO.CPVOlD) &tefnp),TRUE); 

EnableMenuI tem(hwndFrame,COMHANDMSG(&insg)->cmd, FALSE); 
return 0; 
case IDH^ACQUIRE: 

au_ou t pu t ( ha , 0x8 1 88 > ; 

U i r^essageBox ( HUNO.DESICTOP , hwndF rame, 

“Now in passthru mode. \nPosit ion Camera and object.\nSelect OK to freeze 

image.*', 

*• I MAG i nc" , 0 , HB_OK I MB_ I CON ASTER I SK I MB^HOVEAB L E ) ; 
au_output ( ha , 0x8000 ) ; 
return 0; 
case IDM_ACQUIRERGB: 

au_output ( ha , 0x818C ) ; 

U i nHessageBox ( KUND.DE SKTOP , hwndF r ame , 

“Now in RGB passthru mode. \nPosit ion Camera and object.\nSelect OK to freeze 

image.", 

“ 1 MAG i ne“ , 0 , NB_OK | MB_I CONASTER I SK | HB_MOVE ABLE ) ; 
au_output ( ha , 0x8002 ) ; 
return 0; 
case IDM_SAVEAS: 
case IDM^OPEN: 

WinOlgBox(HWND_DESKTOP,hwnd,OpenDlgProc,(KHOOULE)NULL, 

IDD_FILE10,(PV0ID) (COMNANDMSG(&msg)'>cmd)); 
return 0; 
case IDH.BUFFERS; 

W i no I gBox< HWND_DESKTOP , h wnd , Buf 0 1 gPr oc , ( HMOOULE > NUL L , 1 0D_BU F , NULL ) ; 
return 0; 


case IDM^ABOUT: 

UinDlgBox(HWND_DESKTOP,hwnd,AboutDlgProc,(HHOOULE)NULL,IDD_ABOUT,NULL); 
return 0; 

case IDH^THRESHOLD: 

WinShowWindow(WinLoadDlg(HUND_DESKTOP,hwnd,ThreshDlgProc, 

( HHODULE >NULL , I DD^SCROL L , NULL ) , TRUE ) ; 

EnableMenuI tem(hwndFrame,COMMANDMSG(&msg)->cmd, FALSE); 
return 0; 
case 1DM_ENLARGE: 

BNOTB ( en I a rge_ma i n( ) ) ; 
return 0; 
case IDH^ZOOM: 

BNOTB< zoom^ma i n( ) ) ; 



return 0; 
case IDM_ENHRED: 

BMOTB(enhred_mainO); 
return 0; 

case IDM^EQUALIZE: 

BMOTB( equa I i ze_ma i n( ) ) ; 
return 0; 
case IDM_AVERAG£: 

BNOTB(averageO); 
return 0; 
case IDHJNVERSE; 

BNOTB(inverseO); 
return 0; 
case IDM^HOR: 
case lOM^VER: 

temp = COMMANDMSGC&msg)'>cmd - IDM^HOR; 

Ui nShouWi ndou(Ui nLoacD I g( HWND^DESKTOP , hwnd, I ntD I gProc , ( HMODULE >NULL , 
IDO_IMTEGRATE,(PVOID) &tefnp),TRUE); 

EnableHenuI temChwndF rame, COMMANDHSGC&msg) '>ciTid, FALSE ) ; 
return 0; 
case IDM^FOURIER; 

BNOTB(fftO); 
return 0; 
case IDM_BAND_8: 

BNOTB(bfft(8)); 
return 0; 
case IDH_BAN0_16: 

BN0TB(bfft<16)); 
return 0; 
case IDH^ROBERTS: 

BNOTB( roberts^ma i n( ) ) ; 
return 0; 
case IDK_S0BEL: 

BNOT B ( sobe I _ma i n( > ) ; 
return 0; 
case IDH.PREUITT: 

BNOTB(prewltt_fnain<)>; 
return 0; 
case IDM^KIRSCH: 

BNOTBOci rsch_main( ) ); 
return 0; 

> 

break; 

case UM.ERASEBACICGROUND: 
return <MRESULT)1; 

> 

return U i nOef U i ndowP roc ( hwnd , msg , mp1 , mp2 ) ; 

> 

MRESULT EXPENTRY AboutDlgProc(HWND hwnd, USHORT msg, MPARAM mpl, MPARAH mp2) 
i 

switch (msg) i 

case UM.COMHAND: 

8witch(C0HHANDMSG(&msg)->cmd) < 
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case DID.OK: 
case D1D_CANCEL: 

Ui nD 1 snii ssO I g ( hwnd , TRUE ) ; 
return 0; 

> 

break; 

> 

return W i nDef D I gProc( hwnd , msg , mp1 , mp2 } ; 

> 

MRESULT EXPENTRY IntDlgProcCHWNO hwnd, USHORT msg, MPARAM mpl, MPARAM mp2} 
{ 

USHORT t; 

FILE *fp; 
float * intar; 
int cpo8,numl; 

static CHAR *titleC2]=<”Horizontal Integration”, “Vertical Integration”); 
static float *ptrC2]=Cinteg_h, integ_v>; 

HWND win_hwnd ^ Ui hU i ndowFromID (hwnd, I DU_BAR); 
switch (msg) { 

case UH_INITDLG: 

t = *(USHORT *) PVOIDFRO«MP(mp2); 

Uin$etUindowText(UinUindowFromID(hwnd, IOT_INT),”0”); 
WinSetWindowText(UinWindowFromID(hwnd,IDEJNT),“512”); 

UinSetWindowText(UinUindowFromlD(hwnd,FID_TITLEBAR),title[t]); 

UinSetUindowPtr(win_hwnd,0,(PVOID)ptrtt] ); 

Wi nSetUi ndowUShort ( w i n^hwnd, 4 , t ) ; 

WinSendDlgItefTMsg(hwnd,lDSJNT,SBH_SETSCROLLBAR, 

MP FROMSHORT ( 0 ) , MP FROM2S HORT < 0 , 5 11 ) ) ; 

Wi nSend) I gl ten#isg< hwnd, IDS^I NT , SBM_SETPOS , 

MPFR0M2SHORT(0,0),NULL); 

return 0; 
case WM_VSCROLL: 

epos = WinSencDlgIteiT«sg(hwnd,IDS_INT,SBM_QUERYPOS,OL,OL); 

switch (SHORT2FROMMP(mp2)> C 
case SB_L I NEOOWN: 

epos * min(511,cpos+1); 
break; 

case SB_LINEUP: 

epos = maxC0,cpos-1); 
break; 

case SB_PAGEDOWN: 

epos « min(511,cpos'»'16); 
break; 

case SB_PAGEUP: 

epos = max<0,cpos*16); 
break; 

case S6_SLI0ERTRACX: 

epos = SHORT1FROMHP(mp2); 
break; 
default: 

return 0; 

> 

Ui nOueryD I gl temShor t ( hwnd, lOE^l NT , &numl , FALSE ) ; 
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if (cpos'^’tiumi > 512) cpos=512-numl; 

ginSendDl9lt«T#isg<hund.lDS.INT,SBM_SETPOS,MPFROM2SHORT(cpos,0),NULL); 

spr i nt f < buf , ”Xd" , epos ) ; 

Wi nSetD I gl temText ( hwnd, I DT_I MT , buf ) ; 

if (UinSendDlglteirHsg(hwnd,IDB_AUTO,BM_OUERYCHECK,0L,0L)==1 ) 

Wi nPostMsg< hwnd, WM_^COMHAND , MPFROMSHORT (0 IO_OK ) , n|>2 ) ; 
return 0; 
case UM^COMMAND: 

switch(COHMANDMSG(&msg)->cmd) C 
case IDB^CANCEL: 

EnableMenuI tem(hwndFrame, IDM_H0R+ 

UinQueryWindowUShort(win_hwnd,4),TRUE); 

Ui nD i smi ssD I g ( h wnd , T RUE ) ; 
return 0; 
case IDB^OX: 

epos * WinSencDlgIteiT#isg(hwnd, 1DS^1NT,SBM_QUERYPOS,OL,OL); 

WinOueryOlgl teniShort(hwnd, I0E_INT ,&nunil , FALSE ); 
if (cpos+nunl > 512) nunl * 512-cpos; 

BNOTB ( i ntegC Wi nQueryWi ndowUShort ( wi n_hwnd, 4 ) , epos , nunl ) ) ; 

Wininval idateRect( win^hwnd, NULL, FALSE); 
return 0; 
case 1DB_SAVE: 

if (UinDlgBox(HWND_DESKTOP,hwnd.OpenDlgProc,<HMODULE)NULL, 

IDO_FILEIO,(PVOID)COMHANDMSG(&msg)->cind)) i /* save file now */ 

if <(fp = fopcn(szFileName,“w»)) «= MULL) i 

WinMessageBox(HWND_DESKTOP,hwndFrame, "Could not open file!", 

" I HAG i nc” , 0 , MB_OK | MB^I CONEXCLAHAT I ON | HB_MOVE ABLE ) ; 
return 0; 

> 

intar * (float*) UinflueryUindowPtr(win_hwnd,0); 
for <t=0;t<256;t++) 

fprintf(fp,"X.6f\n«,intar(t]); 

fclose(fp); 

> 

return 0; 

> 

break; 

> 

return WinDefDlgProc(hwnd,msg,mp1 ,mp2); 

> 

MRESULT EXPENTRY HstDlgProcCHWNO hwnd, USHORT msg, MPARAM mpl, MPARAM mp2) 
i 

static CHAR *titlct3]=C" Intensity Histogram*', "Saturation Histogram", 

“Hue Histogram”>; 

static float *ptr [3]=Lhisti ,hists,histh>; 

USHORT t; 

HWND win_hwnd; 
switch (msg) { 

case UM^INITDLG: 

win^hwnd = UinUindowFromIO(hwnd, IDW_BAR); 
t » *(USH0RT *) PV0IDFR0MHP(mp2); 

W i nSe t W i ndowT ex t ( W i nU i ndowF romi D ( h wnd , F 1 D^T I T LEBAR ) , t i t I e C t ] ) ; 
WinSetWindowPtr(win_hwnd,0,(PVOID)ptr [tl ); 
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U 1 nSet Wi ndowUShor t ( wi n_hwnd, 4 , t ) ; 

Ui nSetUi rKJowUShort ( wi n^hvind , 6, 8) ; 

W 1 nSe t W i fxloWUShor t ( w i n_h wnd , 8 , 256 ) ; 
return 0; 
case UH.COMTROL: 

win_hwnd = WiniUindowFromID(hwnd,IDW^BAR); 
WinSetUindoviUShort(win_hwnd,6,SH0RTlFR0MMP0np1) -1DBJ6 ♦ 4>; 
return 0; 
case UM^COHMAND: 

suitch(COMHANDMSG(&msg)->cind) i 
case DID.CANCEL: 

uin hwnd * WinUindoHFromID(hwnd,IDW_BAR); 

EnableMenuI teni(hwndFrame, IDM_INT+ 
WinQueryUindowUShort(win_hwnd,4),TRUE); 

Ui nO i smi ssO I g( hund , TRUE ) ; 
return 0; 
case DID.OK: 

win hwnd = UinUindowFromID(hwnd,IDW_^BAR); 
t = WinQueryWindofcAJShort(win_hwnd,6); 
WinSetWindowUShort(win_hwnd,8, 1 « t); 

BMOTB(calc_hist( (float *)WinQu€ryWindowPtr(win_hwnd,0), 
WinQueryWindowUShort(win_hwnd,4),t)); 

WinInvalidateRect(win_hwnd, MULL, FALSE); 
return 0; 

> 

break; 

> 

return WinDefDlgProc(hwnd,msg,mp1,fnp2); 

> 

MRESULT EXPENTRY BufDlgProc(HWND hwnd, USHORT msg, MPARAM mpl, MPARAH mp2) 
int i; 

switch (msg) { 

case UM^INITDLG: 

for (1*0; i <4; i^**') 

if (act^buf & (1«i)) 

WinSendDlgItenWsg(hwnd,lDCJNT>i,BM_SETCHECX.MPFROMSHORT(1),OL); 
return 0; 
case UH.COMMAND: 

8witch(C0MMANDMSG(&msg)*>cmd) { 
case D1D_0K: 

act_buf = 0; 
for (i=0;i<4; i-*'+) 

act_buf 1= (WinSencDlgItenWsg(hwnd,IDCJNT^i,BM_OUERYCHECK, 
MULL, NULL) ? (1«i) : 0); 

if (UinSendDlgItefTMsg(hwnd,lDC^CLR,BM_QUERYCHECX, MULL, MULL)) 
for (i*0; i<4; i-M*) 

if (Uact^buf & (1«i))) au_buf_clear(ha,i); 
UinOismissDlg(hwnd,TRUE); 
return 0; 

> 

break; 

> 
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return WinDefDlgProc<hwnd,msg,mp1,mp2); 

> 

MRESULT EXPENTRY ThreshDlgProc(HWMO hwnd, USHORT msg, MPARAM mpl, MPARAM np2) 
( 

switch (msg) ( 

case UM_INITDLG: 

W1nSendDlgIten«sg(hwnd,IDD_SCROLLBAR,SBM_SETSCROLLBAR, 

MPFROMSHORK thresh), MPFROM2SHORT(0, 255)); 

U i nSendD I g 1 1 enWsg ( h wnd , I DD_SCROL LBAR , SBM_SETPOS , 

MPFROM2SHORT ( th resh , 0 ) , NULL ) ; 
spr i nt f ( buf , , th resh ) ; 

WinSetDlgItemText(hwnd,IDD_SCROLLNUM,buf); 

BNOTB( ini t_th resh ( (unsigned char) thresh)); 
return 0; 
case WM^VSCROLL: 

switch (SH0RT2FR0MMP(mp2)) < 
case SB.LINEDOUN: 

thresh * min(255,thresh+1 ); 
break; 

case SB^LINEUP: 

thresh = max(0,thresh*1); 
break; 

case SB_PAGEDOWN: 

thresh = min(255,thresh<»‘16); 
break; 

case SB_PAGEUP: 

thresh = max(0, thresh- 16); 
break; 

case SB_SLIDERTRAC»C: 

thresh = SH0RT1FR0MMP(mp2); 
break; 
default: 

return 0; 

> 

U i nSendD lgItefnMsg(hwnd,I OD^SCROL LBAR , SBM_SETPOS , 

MPFROM2SHORT ( th resh ,0),NULL); 
spr i n t f ( buf , ”Xd" , t h resh ) ; 

WinSetDlgItemText(hwnd,IDD_SCROLLMUM,buf); 
return 0; 
case UH^COMHAND: 

switch(COMMANDMSG(&msg)->cmd) ( 
case DIO_OK: 

BNOTB(end_thresh(D); 

WinDismissOlg(hwnd,TRUE); 

EnableHenuItem(hwndFrame, IDM_THRESHOLO,TRUE); 
return 0; 
case DID^CANCEL: 

BNOTB ( end^th resh ( 0 ) ) ; 

WinDismissOlg(hwnd,TRUE); 
EnableMenuItem(hwtxiFrame,IDM_THRESHOLD,TRUE); 
return 0; 
case IDB^UPOATE: 

BNOTB(call_thresh( (unsigned char) thresh)); 



return 0; 

> 

break; 

> 

return W i nOef 0 L gProcChund , msg , mp1 , mp2 ) ; 
> 


void far _loadds PMerr(str) 
char *str; 

W 1 nA I arm( HWND_DES»CTOP , WA^NOTE ) ; 

WinMessageBox(HWND_DESKTOP,hwndFrame,str,NULL,0,MB_OK | MB_ICONEXCLAMATIO*l) 

> 

VOID Fi UDirListBox (HUND hwnd, CHAR *pcCurrentPath) 
i 

static CHAR szDrive [] * ” :« ; 

FILEFINDBUF findbuf ; 

HDIR hDir = 1 ; 

SHORT sOrive ; 

USHORT usDriveNum, usCurPathLen, usSearchCount - 1 ; 

ULONG ulDriveHap ; 

DosQCurOisk (AusDriveMum, AulDriveMap) ; 
pcCurrentPath [0] « (CHAR) usDriveNum * (CHAR) 'S' ; 
pcCurrentPath [1] * •:* ; 
pcCurrentPath [2J = *\\* ; 
usCurPathLen » 64 ; 

DosQCurDir (0, pcCurrentPath > 3, &usCurPathLen) ; 

UinSetDlg] temText (hwnd, I0D_PATH, pcCurrentPath) ; 

UinSencOlgIteirMsg (hwnd, IDD.DIRLIST, LH^DELETEALl, NULL, NULL) ; 

for (sDrive * 0 ; sOrive < 26 ; sOrive^-*-) 
if (ulDriveMap & IL « sOrive) 

< 

szOrive [1] = (CHAR) sDrive (CHAR)*A* ; 

UinSendDlgltenMsg (hwnd, IDD.DIRLIST, LM^INSERTITEM, 
MPFROM2SHORT (LIT^ENO, 0), 

MPFROMP (szDrive)) ; 

> 

DosFindFirst (”*.*»*, &hDir, 0x0017, Afindbuf, sizeof findbuf, 
AusSearchCount, OL) ; 

while (usSearchCount) 
i 

if (findbuf. attrFi le & 0x0010 &A 

(f indbuf .achName [0] || f indbuf .achName [1])) 

UinSendDlgltenMsg (hwnd, IDD^DIRLIST, LMJNSERTITEM, 

MPFR0M2SH0RT (LIT^SORTASCENOING, 0), 



HPFROMP (findbuf.achName)} ; 


DosFindNext (hOir, &findbuf, sizeof findbuf, &usSearchCount) ; 

> 

> 

VOID FiUFUeLlstBox (HWND hund) 
i 

FILEFINDBUF findbuf ; 

HDIR hDir s 1 ; 

USHORT usSearchCount > 1 ; 

UinSendDlglteirMsg (hwnd, IDD.FILELIST, LM_DELETEALl, NUIL, NULL) ; 

DosFindFirst &hOir, 0x0007, ftfindbuf, sizeof findbuf, 

&usSearchCount, OL) ; 

while (usSearchCount) 
i 

UinSendDlglteirHsg (hwnd, IDD_FILELIST, LHJNSERTITEH, 

MPFROM2SHORT (LIT_SORTASCENDING, 0), 

MPFROHP (findbuf.achName)) ; 

DosFindNext (hDir, Afindbuf, sizeof findbuf, AusSearch Count) ; 

> 

> 

HRESULT EXPENTRY OpenDlgProc (HUND hwnd, USHORT msg, HPARAM mp1 , MPARAM mp2) 
i 

static CHAR szCurrentPath [80], szBuffer [80] ; 

SHORT sSelect ; 

static (far pascal *sfunc)(int,char *); 
switch (msg) 
i 

case UM^INITDLG: 

FiUDirListBox (hwnd, szCurrentPath) ; 

Fi llFi leListBox (hwnd) ; 

8witch( (int) mp2) i 
case IDH^OPEN: 

sfunc s au_load; 

WinSetDlgl temText(hwnd, IDD_^HEAD,”Open Fi le"); 

U i nSetD I g 1 1 emT ext ( hwnd, D I D_0X, ••Open” ) ; 
break; 

case IDH_SAVEAS: 
sfunc = au^save; 

U i nSetD I g 1 1 emTex t ( hwnd, I DD^HEAD , ”Save F i I e'* ) ; 

U i nSetD I g 1 1 emTex t ( hwnd, D I D_OX , ••Save** ) ; 
break; 

case ]DB_SAVE: 
sfunc s NULL; 

UinSetD Igl tern! ext ( hwnd, IDD_HEAD,**Save to Neural Net File**) 

WinSetDlgl temText(hwnd,DID_OIC, ••Save**); 

break; 

> 

WinSendDlgltemMsg (hwnd, IDD^FILEEDIT, EM^SETTEXTLIMIT, 
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MPFR0H2SH0RT <80, 0), NULL) ; 

return 0 ; 
case UH^CONTROL: 

if (SHORTIFROWP (fflpi) == IDD_DIRLIST || 

SHORT1FROWP (mpl) == IDD^FILELIST) 
i 

sSelect s (USHORT) UinSendDlgltenMsg (hwnd, 

SH0RT1FR0KMP (mpl), 
LM_QUERYSELECTION, OL, OL) ; 

WinSendDlgltemHsg (hvind, SHORT1FROMHP (mpl), 
LM^^OUERYITEMTEXT, 

MPFR0M2SH0RT (sSelect, sizeof szBuffer), 
MPFROMP (szBuffer)) ; 


switch (SHORT1FROMMP (mpD) // Control ID 

< 

case IDD^DIRLIST; 

switch (SHORT2FROMMP (mpD) // notification code 
i 

case LN_ENTER: 

if (szBuffer [0] «* ' •) 

DosSelectDisk (szBuffer [1] * '3') ; 

else 

DosChDir (szBuffer, OL) ; 

FiUDirListBox (hwnd, szCurrentPath) ; 

FillFi leListBox (hwixl) ; 

UinSetOlgl temText (hwnd, IDD^FILEEDIT, *•*') ; 
return 0 ; 

> 

break ; 

case IDD^FILELIST: 

switch (SHORT2FROMMP (mpD) // notification code 
i 

case LN_SELECT: 

WinSetDlgl temText (hwnd, IDD.FILEEDIT, 
szBuffer) ; 

return 0 ; 
case LNJNTER: 

ParseFi leName (szFi leflame, szBuffer) ; 
WinOismissOlg (hwnd, TRUE) ; 
if (sfunc 1- NULL) { 

BNOTB((*sfunc)(ha,szFi leName)); 

> 

return 0 ; 

> 

break ; 


> 



break 


case UM^COMMAND: 

switch (C0MHANDMSG(&fnS9)->cmd) 

< 

case OlOJX: 

WinQueryOlgltemText (hwnd, IDD_FILEEDIT, 

sizeof szBuffer, szBuffer) ; 
switch (ParseFl leName (szCurrentPath, szBuffer)) 
i 

case 0: 

WinAlarm <HWMD_DESKTOP, WA^ERROR) ; 

Fi I IDirListBox (hwnd, szCurrentPath) ; 

Fi UFi leListBox (hwnd) ; 
return 0 ; 

case 1: 

FillDirListBox (hwnd, szCurrentPath) ; 

Fi 1 1 Fi leListBox (hwnd) ; 

UinSetDlgltemText (hwnd, IDD^FILEEDIT, 
return 0 ; 

case 2: 

strcpy (szFileWame, szCurrentPath) ; 
UinOismissDlg (hwnd, TRUE) ; 
if (sfunc Is NULL) < 

BNOTB((*sfunc)(ha,szFi leName)); 

> 

return 0 ; 

> 

break ; 

case DID.CANCEL: 

UinDismissDlg (hwnd, FALSE) ; 
return 0 ; 

> 

break ; 

> 

return WinDefDlgProc (hwnd, msg, mpl, mp2) ; 

> 

SHORT ParseFi leName (CHAR *pcOut, CHAR *pcln) 

< 

/* 

Input: pcOut -- Pointer to parsed file specification, 

pcin -- Pointer to raw file specification. 

Returns: 0 -- pc In had invalid drive or directory. 

1 ** pcIn was einpty or had no filename. 

2 ** pcOut points to drive, full dir, and file name. 

Changes current drive and directory per pcin string. 

V 



CHAR *pcLastSlash, *pcFileOnly ; 

ULONG uIDriveMap ; 

USHORT usDriveNum, usDirLen = 64 ; 

strupr (pcin) ; 

// If input strino is enpty, return 1 


if (pcIn CO] == *\0») 
return 1 ; 

// Get drive from input string or current drive 

if (pcin Cl] *= *:*) 

< 

if (DosSelectDisk (pcin CO] - *5)*)> 
return 0 ; 

pcin 2 ; 

> 

DosQCurOisk (&usDriveNiiii, &ulDriveMap) ; 

♦pcOum = (CHAR) usOriveNum + (CHAR) *3* ; 

*pcOut-^-^ « *:* ; 

*pcOut^+ s *\\» ; 

// If rest of string is empty, return 1 

if (pcin CO] s= *\0*) 
return 1 ; 

// Search for last backslash. If none, could be directory. 

if (MULL « (pcLastSlash = strrchr (pcin, *\\*))) 
i 

if (iDosChOir (pcin, OL)) 
return 1 ; 

// Otherwise, get current dir & attach input filename 

DosQCurDir (0, pcOut, AusDirLen) ; 

if (strlen (pcin) > 12) 
return 0 ; 

if (*(pcOut ♦ strlen (pcOut) - 1) 1= *\\*) 
strcat (pcOut^+, “W”) ; 

street (pcOut, pcin) ; 
return 2 ; 

> 

// If the only backslash is at beginning, change to root 


if (pcin » pcLastSlash) 
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{ 

DosChDir ("\V\ OL) ; 

if (pcin [1] « *\0*> 
return 1 ; 

strcpv (pcOut, pcIn ♦ 1) ; 
return 2 ; 

> 

// Attefnpt to change directory -- Get current dir if OK 
♦pcLastSlash = *\0* ; 

if (DosChDir (pcin, OD) 
return 0 ; 

DosQCurDir (0, pcOut, &usDirLen) ; 

// Append input filename, if any 

pcFileOnly = pcLastSlash 1 ; 

if (*pcFileOnly == *\0*) 
return 1 ; 

if (strlen (pcFileOnly) > 12) 
return 0 ; 

if (*(pcOut ♦ strlen (pcOut) - 1) 1= 'W*) 
strcat (pcOut++, "W") ; 

strcat (pcOut, pcFileOnly) ; 
return 2 ; 

> 

void EnableMenuItem(HWND hwndParent, SHORT sMenuItem, BOOL f Enable) 

< 

HWND hwncBienu = WinWindowFromID(hwndParent,FID_MENU); 

U i nSencHsg ( h wncWenu , MM^SET I TEHATTR , MP FROM2SHORT ( sMenu 1 1 em, TRUE ) , 
MPFROM2SHORT(MIA_DISABLEO, fEnable ? 0 : MIA^DISABLED)); 

> 


no 


/* me INTEG.C V 

#include <stdio.h> 

^include ”util.h“ 

#define HEIGHT 512 
#define WIDTH 512 

#define max2(8,b> ((a) > (b)) ? (a) : (b) 

extern unsigned char huge ipIcIHEIGHT] [WIDTH] ; 

extern unsigned char huge opic [HEIGHT] [WIDTH] ; 

float 1 nteg_h [2561 , 1 nt eg_v [256] ; 

void IntegvC unsigned int, unsigned int); 

void InteghC unsigned int, unsigned Int); 

void 1nteg(unsigned int, unsigned Int, unsigned int); 

extern int ha; 

void integ(or,start,num) 
unsigned Int or, start, nun; 
i 

switch (or) i 
case 0: 

1 ntegh( start , nun) ; 
break; 
case 1: 

i ntegvC star t , nun) ; 
break; 

> 

> 

void InteghC unsigned int start, unsigned int nun) 

unsigned int i,j; 

float mx= (float) 0.; 

au_readscr(ha, ipic,0); 

for (i=0;i<256;i+-«') integ_h[i]= (float) 0.; 

for (i«0; i<512;i++) 

for ( j=start; j<start+nun; j++) integ_h[1/2]^= (float) ipic[i][j]; 
for (i-0;i<256;i++) mx = max2(mx, integ_h[1] ); 
if (mx (float) 0.) 

for (i=0;i<256;i*+) 1nteg_h[i] /= mx; 

> 

void integvC unsigned Int start, unsigned int nun) 

< 

unsigned int i,j; 
float mx= (float) 0.; 
au_readscr(ha, ipic,0); 

for (i=0;i<256;i^+) integ_v[i)« (float) 0.; 
for (i=start;i<start+nun;i++) 
for ( j=0; j<512; ]♦♦) 

integ_v[j/2]*= (float) ipic[i][j]; 
for (i=0;1<256;H+) mx * max2(mx, integ_v[i] ); 
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if (mx 1= <float)0.) 

for (i=0;i<256;i++) integ_v[i] /= mx; 

> 
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/* File MAPC.C */ 

; Static Name Aliases 

! 

TITLE mapn.c 

NAME mapn 

.286p 

.287 

MAPNJEXT SEGMENT WORD PUBLIC ’CODE* 

MAPN_TEXT ENDS 

_DATA SEGMENT WORD PUBLIC 'DATA' 

_0ATA ENOS 

CONST SEGMENT WORD PUBLIC 'CONST* 

CONST ENDS 

_BSS SEGMENT WORD PUBLIC 'BSS' 

_BSS ENDS 

DGROUP GROUP CONST, _BSS, _DATA 

ASSUME CS: MAPN^TEXT, DS; DGROUP, SS: DGROUP 
EXTRN _acrtused:ABS 
EXTRN _AHINCR:FAR 
MAPN_TEXT SEGMENT 

ASSUME CS: MAPN.TEXT 

; Line 4 

PUBL I C _mapn 

_mapn PROC FAR 

enter WORD PTR 22,0 
push di 

push si 

; Line 6 
; pic = 6 

; map = 10 

; li s -8 

; i = -10 


tpic = 
10 

-4 


mov 

WORD PTR [bp-10],0 

;i 

mov 

WC*0 PTR [bp-81,0 


mov 

WORD PTR [bp-61,1 


mov 

ax,WORD PTR [bp«^61 

;pic 

mov 

dx,WORD PTR [bp+81 


mov 

WORD PTR [bp-201, ax 


mov 

WORD PTR [bp-181,dx 


mov 

di,4 


add 

WORD PTR [bp- 101, 4 



SL20002: 

; Line 11 

mov ax, WORD PTR [bp-20] 

mov dx,UORD PTR [bp-18] 

mov WORD PTR [bp-4], ax ;tpic 

mov WORD PTR [bp-2] ,dx 

; Line 12 


mov 


WORD PTR [bp-14],0 
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mov 

WORD PTR [bp-12],1 


mov 

WORD PTR [bp-16],ds 


les 

si, DWORD PTR [bp-4] 

;tpic 

; Ids 

si, DWORD PTR [bp-4] 

;tpic 

Ids 

bx, DWORD PTR [bpflO] 

;map 

fflOV 

cx,65535 


$120000: 

; Line 13 

mov 

at, BYTE PTR es;[si] 

;al = *tpic 

xlatb 


;al 3 maplal] 

mov 

BYTE PTR as: [si], at 

;*tpic = al 

inc 

si 


loop 

SL20000 


mov 

al,BYTE PTR es:[si] 


xlatb 

mov 

BYTE PTR es:[si],al 


mov 

WORD PTR [bp-4], si 

;tpic 

mov 

WORD PTR [bp-2], as 


mov 

ds,WORD PTR [bp- 16] 


ASSUME 

DS: NOTHING 


; Line 14 

mov 

ax,OFFSET _AHINCR 


add 

WORD PTR Cbp'ISl.ax 


dec 

di 


jnc 

SL20002 


; Line 15 

pop 

si 


pop 

di 


leave 

ret 

_mapn ENDP 

MAPN^TEXT 

ENDS 


END 



File PROCESS. C 






#include «util.h“ 

#define LORED 255 
#define HIRED 240 
#define SATJHRESH 32 

unsigned char huge tpicC512] [512] ; 

unsigned char huge ipic[512] [512] ; 

unsigned char huge opic[512] [512]; 

static unsigned int ihist[256]; 

static unsigned char map [256]; 

void zooffl_main(void); 

void enlarge_main(void); 

void enhred_main(void); 

void equal ize_main( void); 

void zoom(int); 

void enlarge(int buf); 

void init_thresh<unsigned char th); 

void cal l_thresh( unsigned char th); 

void end_thresh(int); 

void average(void); 

void mapnCunsigned char huge unsigned char *); 

void copyCunsigned char huge unsigned char huge *); 

extern int ha; 

extern unsigned int act_buf; 

void zoom_main(> 
i 

int i; 

/* Initialization */ 
for (i=0;i<4;i++) 

if (act_buf & (1«i)) 
zoom(i ); 

> 


void zoom(buff) 
int buff; 
i 

int ij,x,y; 

au_r eadsc r ( ha , i pi c , buf f ) ; 

for (i»0 , x«128 ;i<512;i^»2 , x-m-) 

for (j=0 , y=128; j<512; j+=2 , y++) { 
opic[i][j] * ipiclxHy]; 
opic[i] [j*1]= (ipic[x] [y]^ipic[x] [y+1])/2; 
opic[i+1] [j]= <ipic[x] [y]-**ipic[x^1] [y])/2; 
opic[i+1] [j^1]= (ipic[x] [y]+ipic[x+1] [y+1] )/2 
> 

au_wr i tesc r ( ha « opi c , buf f ) ; 

> 


void enlarge^mainO 
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{ 

int i; 

!* Initialization V 


for (1*0; 1 <4; i-*^) 

If (act^buf & (1«D) 
enlarged); 

> 

void enlarge(buf) 
int buf; 

< 

int ij; 

unsigned char huge *topic1 * (unsigned char huge *) opicIO); 
unsigned char huge *topic2 » (unsigned char huge *) opicH); 
unsigned char huge *tipic * (unsigned char huge *) &ipic[128] [128] ; 
au_readscr ( ha , i pi c , buf ) ; 

for ( i =0; i <256; i , topi c1 =topi c2 , topi c2+=51 2 , t i pi c^=256) 
for ( j*0; j<256; j^>,topic1+=2,topic2+=2) 

♦(topici) * *(topic1+1) * *(topic2) * *(topic2+1) * *(tipic+^); 
au_wr i t esc r ( ha , opi c , buf ) ; 

> 

void enhred_main() 

{ 

long li; 

unsigned char huge *topic » (unsigned char huge *) opic; 

unsigned char huge *tipic * (unsigned char huge *) ipic; 

unsigned char huge *ttpic * (unsigned char huge *) tpic; 

au^buf _c I ear ( ha , 3 ) ; 

au_readscr(ha,tpic,0); 

au_readscr ( ha , i pi c , 2 ) ; 

au_readscr (ha , opi c , 1 ) ; 

for (li«0L;li<262144L;li++,tipic'M*,topic++,ttpic+^) 

if (((*tipic <= LORED) || (*tipic >= HIRED)) (*topic > SAWHRESH)) { 
*ttpic * 255; 

*tipic * 255; 

•topic * 255; 

> 

else i 

•ttpic = 0; 

> 

ttpic * (unsigned char huge •) &tpic[480] [0]; 
for ( I i *0L ; I i <32L*5 1 2L ; I i , 1 1 pi C++ ) 

•ttpic * 0; 

au_wri tescr(ha,tpic,0); 
au_writescr(ha,ipic,2); 
au_uritescr(ha,opic,1); 

> 

void equalize_inain() 

< 

const unsigned int bins*256; 
unsigned int binpos=0; 


unsigned int i; 
unsigned long sliipOL; 
unsigned long maxbinva 1^10241; 
unsigned long li; 

const unsigned long binwidth»1024L; 
unsigned char huge *tp « ipicCO]; 

r 

maxbinvalsbinwidth-(512l*512L)/(unsigned long) bins; 

au_readscr(ha,1pic,0); 

for (i*0;i<256;i-M») ihist[i]>0; 

for (li*0L;li<262144L;li‘*^,tp*^') 

‘M'ihi St [(unsigned int)(*tp) ]; 
for (i*0;i<256;i-M-) < 

while ((sunH^C unsigned long)ihist [i] ) > maxbinval) f 

binpos-M-; 

maxbinval binwidth; 

> 

map[i]s (unsigned char) (binpos*256/bins); 
sum (unsigned long) ihist[i]; 

> 

tp s ipicCOl; 
mapn(tp,map); 
au_wr i tcsc r ( ha , i pi c , 0 ) ; 

} 


void init^thresh(th) 
unsigned char th; 

< 

unsigned int i; 
unsigned char map [256]; 

unsigned char huge *tipic = (unsigned char huge *) tpic; 
unsigned char huge *topic * (unsigned char huge ♦) opic; 
au_readscr(ha,tpic,0); 
copyC t i pi c , t opi c ) ; 

for (i=0;i<256;i-*-+) mapCi] = (unsigned char) ((i<th) ? 0:255) 

mapnC topic, map); 

au_writescr(ha,opic,0); 

> 

void cal l_thresh(th) 
unsigned char th; 
i 

unsigned int i; 
unsigned char map [256]; 

unsigned char huge *tipic * (unsigned char huge *) tpic; 
unsigned char huge *topic * (unsigned char huge *) opic; 
copy( t i pi c , topi c ) ; 

for (i=0;i<256;i*M-) map[i] * (unsigned char) ((i<th) 7 0:255) 

mapn( topic, map); 

au_writescr(ha,opic,0); 

> 



void end_thresh(f lag) 
int flag; 
i 

if (Iflag) 

8u_wr i tescr(ha, tpi c, 0) ; 

) 

void average(void) 
int ij,t; 

au_readscr(ha, ipic,0); 
for (is1;i<511;i*M») 

for (j=1;j<511;j++) { 

t * (lnt)ipic[i-1] [j-l]+(int)ipicCi-1) [j)+(int)ipicli] Ij-1] ; 
t (int)ipic[i) [j]+(int)ipic[il Cj+1]+(int)ipic[i+13 [j3+(int)ipic[i+1] [j+1] 
opictilCj] = (unsigned char) (t/7); 

> 

au_writescr(ha,opic,0); 

> 

void inverse(void) 
i 

int i,j; 

au_readscr(ha, ipic,0); 
for (i*1; i<511;i^^) 

for (jsl; j<511;j+-»') i 

opic[i][j] = (unsigned char) (255-ipictiHj] ); 

> 

au_writescr(ha,opic,0); 

> 
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/♦ File TEXTURE. C V 

#include "util.h” 

^include <m8th.h> 

#define max2(a,b) ((a) ><b)) ? (a) : (b); 

#define Height 512 
#define Width 512 

extern unsigned char huge ipic [Height] [Width] ; 

unsigned char huge Ak[Height] [Width]; 

unsigned char huge lk[Height] [Width]; 

float hist [256] , h i st i [256] , h i s th [256] , h i s ts [256] ; 

extern int ha; 

void one_iter(void); 

void threshold_regions(void); 

void averagejDicCvoid); 

unsigned char ret_region_aver( unsigned int, unsigned int, unsigned int); 

unsigned int pi ck_region< unsigned int, unsigned int); 

unsigned int pi ck_spec_region( unsigned int, unsigned int); 

unsigned char avert unsigned int, unsigned int); 

unsigned int homogent unsigned int, unsigned int); 

void calc_hist(f loat *, unsigned int, unsigned int); 

void main(void); 

void doitall(void); 

#ifdef DSGSDFGDFG 
void doi tall (void) 
i 

one_iter(); 
one_iter(); 
th resho I d_reg i ons ( ) ; 

> 

void one_iter(void) 

C 

unsigned int i,j; 
au_readscr(ha,ipic,0); 

8veragej3ic(); 

for (i=0; i<Height; i^*»>) 

for(j=0; j<Width; j++) 

Ik[i][j] s (unsigned char) 0; 

for (i=4; i<Height-4; i++) 

for(js4; j<Width*4; ]♦+) 

Ik[i][j] * ret_region_aver(i, j, picker eg ion(i,j)); 
for (i*0; i<Height; i++) 
for(JsO; j<4;j-*-^) C 

Ik[i][Jl = ret_region_aver(i,i,pick_spec_region(i, j)); 

lk[i] [Width-1- j] = ret^region aver(i, Width-1- j, pick spec region(i , Width- 1- 

> 

for (i=0; i<4;i^+) 
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for (j=4; j<Width-4; j++) ( 

1 k [ n [ j ] * ret^reg i on_avcr( i , j , pi ck_spec_reg i on ( i , j ) ) ; 

Ik[Height-1-i] [j] = ret_region_aver(Height-1-f , j,p^ck_spec_^eg^on(Height-1- 

^j»; 

> 

au_writescr(ha,lk,0); 

> 

#endif 

void calc_hist(hist,buf,hist_level) 
float *hist; 

unsigned int buf,hist_level; 

< 

float mx»(f loat>0.; 

unsigned char huge *tp * (unsigned char huge *) ipic; 
int i; 

unsigned long li; 
au_readscr(ha,ipic,buf); 
for (i«0;i<256;i^^) hist[i] ■ (float)O.; 
for (li=0;li<262U4L;li-M>,tp*-+) 

hist[ ((unsigned int) (*tp) » (8-hist_level))l ♦= (float) 1.; 
for (i=0;i<256;i-*-+) mx = max2(mx,hist[i]); 
for (i=0;i<256;i++) hist[i) /= mx; 

> 


#ifdef DFGOFGOFGFDGFDG 
unsigned char rct_region_aver(i , j,d) 
unsigned int i,j,d; 
i 

switch (d) < 

case 0: 

return(aver(i, j)); 
break; 

case 1: 

return(aver( i+2, j+2)); 
break; 

case 2: 

return(aver(i-2, j+2)); 
break; 

case 3: 

return(aver(i-2, j-2)); 
break; 

case 4: 

return(aver(i+2 J-2)); 
break; 

> 

> 

unsigned char aver(x,y) 
unsigned int x,y; 
i 

unsigned int t; 

t = ipic[x-2] Cy-2)+ipicCx-1) Cy-2)+ipiclx] (y-2)+ipic[x+1) Cy-2)+ipicCx+2) Cy-2); 
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t ♦* ipictx-23 Cy-13-M‘pic[x-1] Cy-13^ipiclx3 Iy-l3+ipic[x+13 Cy-13-M‘picCx+23 Cyl]; 
t ipicCx-23 Cy3+ipicCx-13 Cy3+ipic£x] ty3+ipic[x+1] Cy]+ipicCx+2] ty3; 
t ♦= ipicCx-23 tv^13*M’picCX“13 £y^l3+ipicCx3 £y+13*M*pic£x+1] £y+l3+ipic£x>23 £y+13; 
t ♦* ipic£x-23 £y+23+ipic£x-13 £y+23+ipic£x3 £y+23+ipic£x+1] £y+2]+ipicEx+2] £y+23; 
t /* 25; 

returnC (unsigned char) t); 

> 


void threshold_regions(void) 


unsigned int i,a,b, jsO.ksO; 

unsigned int valleys £103, peaks £103; 

unsigned int hist_level * 8; 

unsigned int nuii_hi8t_level8= 1 « hist_level; 

unsigned int level, I ow_ptBO,high_ptsO, modes 1; 

unsigned char map £2563; 

float lou_vals (float) 1., high_val= (float) 0.; 

calc_hist(hist,0,hist_level); 

au^readscr ( ha , i pi c , 0 ) ; 

vallcys£j-M-3 = 0; 

for (i=0; i<n(jn_hist_ levels; i++) i 
suitch(mode) i 

case 0: /* looking for minimum V 
if (hist £i3<=low_val) { 

low^val = hist£i3; 
lowjDt * i; 

> 


else if (hist(i3 > (float)3.*low_val) { 
valleysEj++3 * low_pt; 
mode s 1; 

high^val * hist£i3; 
high^pt s i; 

y 

break ; 

case 1: /* looking for maximum V 

if (histEi] >s high^val) i 
high^val * hist£i3; 
highjDt s i; 

) 

else if (hist£i3*3, < high^val) < 
peaks Ek-*'-*-) * highjjt; 
mode s 0; 
low^val s hist£i3; 
low_pt s i; 

> 

break; 

) 

> 

switch (mode) < 
case 0: 

valleys£j++3 = i-1; 

break ; 

case 1: 


peaks £k++3 = high_pt; 
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valleys [j++] = i-1; 
break; 

> 

b s 0; 

for (i*1;i<j;i++) < 

level = peaks Ci- 13 *256/ rHJii_hist_levels; 
for (a=b; a<(valley8[i]+1)*256/nun_h1st_levels; a+-**) 
mapta] « (unsigned char) level; 

b * a; 

> 

for(i»0;i<Height 

for( j=0; j<Width; j-M-) { 

k « (int) (ipicCi3[j]); 
ipicEi] (j) s map[k3; 

> 

au_writescr(ha, ipic^O); 

> 


unsigned int pick_spec_region(x,y) 
unsigned int x,y; 

< 

unsigned int mnh,th; 
unsigned int dir; 
wnh = 65535; 

if ((x>*2) && (x<*253) && (y>s2) && (y<=253)) C 
if ((th=homogen(x,y)}< mnh) < 
nnh * th; 
dir = 0; 

> 

> 

if ((x<=251) && (y<=25D) { 

if ((th=homogen(x+2,y+2))< mnh) < 
mnh = th; 
dir = 1; 

> 

> 

if ((x>=4) && (y<=25D) { 

if ((th=hoiTK>gen(x-2,y+2))< mnh) < 
mnh = th; 
dir = 2; 

> 

> 

if ((x>s4) && (y>«4)) < 

if ((thshonK>gen(x-2,y-2))< mnh) i 
mnh s th; 
dir = 3; 

> 

> 

if (<x<=251) && (y>=4)) < 

if ((th=homogen(x+2,y-2))< mnh) < 
mnh = th; 
dir = 4; 

> 
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> 

return(dir); 

> 

unsigned int pick_rcgion(x,y) 
unsigned int x^y; 
i 

unsigned int mnh,dir,th; 
mnh ■ honx>gen(x,y); dir * 0; 
if <(th=homogen<x+2,y+2))< wnh) < 
imh = th? 
dir = 1; 

> 

if ((th=homogen(x-2,y+’2))< mnh) i 
mnh » th; 
dir = 2; 

y 

if ((th=homogen(x-2,y-2))< mnh) { 
mnh = th; 
dir = 3; 

> 

if ((th=homogen(x+2,y*2))< mnh) C 
mnh * th; 
dir = 4; 

> 

return(dir); 

> 

unsigned int homogen(x,y) 
unsigned int x,y; 
i 

unsigned int xp = x*1 ,xro=x-l,yp=y+1,ym=y-1; 

return(abs( (int) <Ak[xp) lypl+Aktxm] Cyp)-Ak[xp) £ym]-Ak(xnO tynO) ) ♦ 
abs( (int) (Ak[xp] [yp]^Ak[xp] [ym] *Ak[xm] [yp]-Ak[xm] [ym] ))); 

> 

void average_pic() 

< 

unsigned int i,j; 

unsigned int te; 

for (i*0; i<Height; i+^) 

for ( j=0; j<Width; j-^+) 

AkCiJCj) = (unsigned char) 0; 
for (i=l; i<Height-1; i++) i 

te * ipicCi-11 C0)+ipic(i) t01+ipic[i+1) [0] ♦ ipicCi-1) t1) + ipic[i) [1)+ipicCi-»-l] ED ; 
te /= 6; 

AkEilEO] s (unsigned char) te; 

te = ipicCi-1) CWidth-1)+ipic[i3 [Width-D + ipicEi+D CWidth-1); 
te ♦= ipicEi-1) EUidth*21-»*ipicEi] EWidth-2]+ipicEi*»‘1] EWidth-2] ; 
te /= 6; 

AkEi] EWidth*1] ^ (unsigned char) te; 
te = ipicEO] Ei-ll^ipicEOHil + ipicEO) Ei*1] ; 
te = ipicED Ei-1]-*-ipicE1) Eil*»*ipicE1] Ei+1]; 
te /= 6; 
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Ak[0][i] s (unsigned char) te; 

tc s ipic(Height-1] [i-1]+ipic(Height-D (iJ^ipicCHeight-1] [i+1] ; 
te * ipic[Height-2] (i-l]+ipicCHeight-23 Cil-^ipic[Height-2] [i+1]; 
te /* 6; 

Ak(Height-1Hi] « (unsigned char) te; 
for( jal; j<Uidth-1; j-M-) { 

te ■ ipic[i-l] Cj-11 + ipicCi-l] [jl^ipicCi-1] Ij+1); 
tc ipic[i] Cj-1]*ipicCi] [jJ+ipicCi) 

te ipic[i+l] [j-1)+ipic[i+1Hj]+ipicCi+1] [j+13; 

te /- 9; 

Ak[i][j] « (unsigned char) te; 

) 

> 

te « (ipiclO](03+ipic[03[13+ipicC1)[0)+ipic(1]C1])/4; 

Ak[0][0] s (unsigned char) te; 

te = (ipictO) (Width- 1)+i pic CO] CUidth-2)+ipic[l] CWidth-l]-*-ipicCl3 [Width-2] )/4; 

AkCO] [Width- 1] = (unsigned char) te; 

te = ( ipic [Heigh t-1] [0]+ipicCHeight-1] [13+ipic [Height-2] CO] +i pic [Height -2] C1])/4; 

Ak [Height- 1] [0] - (unsigned char) te; 

te = (ipic[Height-1] [Width-1]*ipic[Height-1] [Width-2]+ipic[Height-2] [Width-1]+ipic[Height- 
2] [Width-2] )/4; 

Ak[Height-1] [Width-1] = (unsigned char) te; 

> 

#endif 



APPENDIX C 


Source Code Listing for Device Drivers 
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# MAKE.MC 

# Makefile for BP 

# 

CC « cl 

#CFLAGS= /AL /FPi87 /G2 /Zi /Od /W3 
CFLAGS* /AL /FPi87 /G2 /Oailt /Gs /W3 

.c.obj; 

S(CC) /c $<CFLACS) $*.c 
.obj.exe: 

# link /CO /NOD $*.obj util.obj,$*.exe,nul.inap, llibc7p c:\os2\doscal Is. lib, def.def 
link /NOO $*.obj utU.obj,$*.exe,nul .map, llibc7p c:\os2\doscalis. lib, def.def 

•c.exe: 

S(CC) /c S(CFLAGS) $*.c 

# link /CO /NOO S*.obj,$*.exe,nul.map, llibc7p c:\os2\doscaUs.lib util. lib,def .def 
link /NOO S*.obj,$*.exe,nul.map,llibc7p c:\os2\doscalls.lib util. lib, def .def 
del S*.obJ 


load.exe: ioad.c 
save.exe: save.c 
enlarge.exe: enlarge. c 
zoom.exe: zoom.c 
equal.exe: equal. c 
enhred.exe: enhred.c 
clrbuf.exe: clrbuf.c 
acquire.exe: acquire.c 
on.exe; on.c 
off.exe: off.c 
paint.exe; paint.c 
dual.exe: dual.c 


half.exe: half.c 



File ACQUIRE. C 


n 




#include <stdio.h> 

^include <conio.h> 

^include "util.h” 
void main(void); 

void mainO 

int ha; 

au_init(&ha); 

au_out put ( ha , 0x81 88) ; 

printf(«Now in passthru mode.\n»); 

prfntf(”Hit any key to freeze image. \n*'); 

getchO; 

au_out put ( ha , 0x8000 ) ; 
au_end(ha); 

> 



File CLRBUF.C 


7 




#include <stdio.h> 

#include <stdlib.h> 

#ir>clude ■•util.h” 
void mainCint.char *□); 

void main(argc,argv} 
int argc; 
char *argv[]; 
i 

int ha; 

int buff^i; 

if (argc«1) exit(O); 

au_init(&ha); 

for (is1;i<argc;i++) < 

buf f=atoi (argvCi] ); 

printf (“Clearing buffer Xd. . ,\n'*,buff ) 
au_buf_clear(ha,buff ); 

> 

au_end(ha); 

> 
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/* File OUAL.C V 

#incliide <stdio.h> 

#include <stdlfb.h> 

#include "util.h" 

typedef int WORD; 
typedef long DWORD; 

unsigned char huge pic[512] C512] ; 
void main( int, char *[]); 

void main(argc,argv) 
int argc; 
char *argv[); 

FILE *fpl,*fp2; 
int i,j,ha; 
if (argc 1*3) exit(O); 
au_init(&ha); 

printf ("Loading Xs. , ,\n",argv(1] ); 
printf( "Loading Xs. . .\n",argv[23 ); 
fp1=fopen(argvIl3,"rb"); 
fp2*fopen(argv[23 /*rb">; 
f seek( f pi , 51 2L , SEEK^SET ) ; 
f seek ( f p2 , 5 1 2L , SEEK^SET ) ; 
for (i*0;i<512;i-M-) 

for (j=0;j<256; j++) i 

pic(i3Cj3 * (unsigned char) fgetc(fpl); 
fgetc(fpl); 

pic[i3 [256^ j 3 -(unsigned char) fgetc(fp2); 
fgetc(fp2); 

> 

au_writescr(ha,pic,0); 
for (i=0;i<512;i^-*-) 

for (j*0; j<256; ]♦♦) < 

pic[i3Cj3 * (unsigned char) fgetc(fpl); 
fgetc(fpl); 

picCi) C256+j)s(unsigned char) fgetc(fp2); 
fgetc(fp2); 

> 

au_writescr(ha,pic,1); 
for (i=0;i<512;i++) 

for (j=0;j<256;j++) i 

pic[i3C]3 = (unsigned char) fgetc(fpl); 
fgetc(fpl); 

picCi) C256+j)s(unsigned char) fgetc(fp2); 
fgetc(fp2); 

> 

au_writescr(ha,pic,2); 

fclose(fpl); 

fclose(fp2); 
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au_end(ha); 

> 
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/* File ENHRED.C V 

#include <stdio.h> 

#1nclude **uti I .h” 

#define LORED 40 
#define HIRED 185 


int status, colorval; 

int values[256], actrgn[41; 

char filename [20], comment [100] , ch; 

unsigned char huge huepic[512] [512]; 

unsigned char huge satpic[512] [512]; 

void main(void); 

void mainO 

int i,j,ha; 

au_init(&ha); 

au_buf _c I ea r ( ha , 3 ) ; 

au_readscr ( ha , huepi c , 2 > ; 

au_r eadsc r ( ha , sa tpi c , 1 ) ; 

printf (“Enhancing red areas of image. - .\n“); 

for (i*0; i <480; i *►•*•) 

for (i=0; j<512; j-M-) < 

if (<huepic[i] [j] <=LORED) || (huepic[i] [j]>=HIRE0)) { 
huepic[i] [j]=255; 

> 

else 

satpicli] [j]*0; 

> 

au_wr i tescr(ha,satpic, 1 ); 
au_writescr(ha,huepic,2); 

printf (“Bright red areas show evidence of liquid penetrant.\n“); 
au_end(ha); 

> 


File ENLARGE. C 


7 


#include <stdio.h> 

#include <conio.h> 

#include "util.h“ 

unsigned char huge ipic[512] [512] ; 
unsigned char huge opic[512] [512]; 
void main(void); 
void enlargeCint buf); 

int ha; 

void mainO 
i 

int i; 

/* Initialization */ 

au_init(&ha); 
for (i=0;i<3;i++) i 

printf(**Enlarging buffer %d...\n*',i) 
enlarge(i); 

> 

au_end(ha); 

> 

void enlarge(buf) 
int buf; 
i 

int ij; 

unsigned char ch; 
au_readscr(ha,ipic,buf }; 
for (i=0;i<512;i^=2) 

for < j=0;i<512;i+=2) i 

ch = ipic[128+i/2] [128+J72]; 
opic[i] [j] = ch; 
opicCi] [j+11 s ch; 
opic[i+1] [j] = ch; 
opicCi+1] [j+1] = ch; 

> 

au_writescr(ha,opic,buf ); 

> 



/* File EQUAL. C V 

#include <stdio.h> 

#include <stdlib.h> 

#include "util.h” 

unsigned char huge ipicC512] C512); 
void mBin(void); 

void mainO i 
int ha; 

unsigned int bin$s256; 
unsigned int hist [256] ,j,map [256] ,binpos, i; 
unsigned long binwidth,sun,inaxbinval,li,lj; 
au_init(&ha); 
au_readscr(ha, ipic,0); 
printf (“Initial izing...\n«); 
for (i=0; i<256;i^^) 
hist Ci)=0; 

maxbinval=binwidth=(512L*512L)/(unsigned long) bins; 

binpossO; 

sum=0L; 

printf (“Entering histogram. . .\n“); 
for (li=0L;li<512L;li-»-^) < 

for (lj=0L;lj<512L;lj++) { 

i=(unsigned int)ipic[li][l j] ; 

♦♦histCi] ; 

> 

> 

printf( “Entering mapping. . ,\n“); 
for ( j=0; j<256; < 
loop: 

if ((su7H-( unsigned loog)hist [j] ) <= maxbinval) C 
map [ J) =bi npos*256/bi ns ; 
sum ♦= (unsigned long) histCj); 

> 

else < 

bi npos>+; 

maxbinval += binwidth; 
goto loop; 

> 

> 

printf (“Entering remapping. . .\n“); 
for (li*0L;li<512L;li+^) { 

for (lj=0L;l j<512L;lj-M.) < 

i = (unsigned int) ipicCli] [1 j) ; 
ipicriiHlj] = (char) mapCi]; 

> 

> 

au_wr i tescr(ha, ipic,0); 
au_end(ha); 

> 
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/* File HALF.C V 

#include <8tdfo.h> 

#include <stdlib.h> 

#include “util.h** 

typedef int WORD; 
typedef long DWORD; 

unsigned char huge pic [512] [512]; 
void main( int, char *[]); 

void main(argc,argv) 
int argc; 
char *argv[]; 
i 

FILE Mp1,*fp2; 
int ij,ha; 
if (argcls3) exit(O); 
au_init(&ha); 

printf (■^Loading Xs. . .\n“,argv[1] ); 
printfC “Loading Xs. . .\n“,argv[2] ); 
fp1=fopen(Brgv[lJ ,••rb“); 
fp2=fopen(argvC2] /‘rb**); 
f seek ( f pi , 5 1 2L , SEEK^SET ) ; 
f seek ( f p2 , 5 1 2L , SEEK^SET ) ; 
for (i*128;i<384;i-M-) < 

for (j=0; j<256; j++) { 

pic[i][j] E (unsigned char) fgetc(fpl); 
fgetc(fpl); 

pic(i) [256+j]*(unsigned char) fgetc(fp2); 
fgctc(fp2); 

> 

f seek ( f pi , 5 1 2L , SEEK_CUR ) ; 
fseek(fp2,512L,SEEK_CUR); 

> 

au_writescr(ha,pic,0); 
for (i*128;i<384;i^+) { 

for (j=0; j<256; j-M*) { 

pic[i][j] £ (unsigned char) fgetc(fpl); 
fgetc(fpl); 

pic[i] [256+j]=(unsigned char) fgetc(fp2); 
fgetc(fp2); 

> 

f seek( f pi , 5 1 2L , SEEK_CUR ) ; 
f seek( f p2. 512L , SEEK_C:UR ) ; 

> 

au_writescr(ha,pic,1); 
for (is128;i<384;i+>) c 

for (j=0;j<256; ]♦+) i 

pic[i][j] = (unsigned char) fgetc(fpl); 
fgetc(fpl); 



picCiJ C256+J)s(unsigned char) fgetc(fp2) 
fgetc(fp2); 

> 

f seek ( f pi , 5 1 2L , SEEK^CUR ) ; 
f seck( f p2 , 5 1 2L , SEEK^CUR ) ; 

> 

au_wri tescr(ha«plc,2); 
fclose(fpl); 
fclose(fp2); 
au_cnd(ha); 

> 
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/♦ File LOAD.C V 

#include <stdio.h> 

#include <stdllb.h> 

#include <string.h> 

#include "util.h" 

typedef int WORD; 

typedef long DWORD; 

void main(int argc, char *argv[]); 

void inain(argc,argv) 
int argc; 
char *argv[}; 

C 

char f name [130]; 
int ha; 

if (argc==2) strcpy(fnamc,argv[1] ); 
else C 

printfC “Enter file to load: "); 
scanf ( “Xs" , f name) ; 

> 

au_init(&ha); 

printfC “Loading Xs...\n”,fname); 
au_ I oad ( ha , f name ) ; 
au_end(ha); 

> 
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/* File OFF.C V 

#include "util.h" 
void inain(void); 

void inainO 

< 

int ha; 
au_init(&ha); 
au_output < ha , 0x0000 ) ; 
au_end(ha); 

> 
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/• File ON.C V 

#include “utiL.h" 
void Riain(void); 

void mainO 
i 

int ha; 
au_init(&ha); 
au_output ( ha , 0x8000 } ; 
au_end(ha); 

> 
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/* File PAINT. C V 

#define INCL^BASE 
#define INCL.DOSDE VICES 
#include <std1o.h> 

#include <stdlib.h> 

#include <os2.h> 

#include “util.h** 

typedef int WORD; 
typedef long OUORD; 

unsigned char huge pic[512][512] ; 

MOUE VENT INFO MouEvent; 

mainO 

i 

int ij; 

char OutStr[40]; 

USHORT Cells0x0720; 

HMOU MouKandle; 

USHORT status; 

SHORT UaitOption^l; 
int pi an [4]; 
au_ini t(); 

pi an [0] spl an M ] «pl an t2] -pi an [3] ^0; 

s t a t us=MouOpen( OL , AMouHand I e ) ; 

au_set_ovl _plns(plan); 

au_se t_ov I _c I rs ( p I an ) ; 

au_se t_cons t ( 3 , 255 ) ; 

getch(>; 

if (status) < 

printf(“MouOpen Failed.\n"); 
exit(l); 

> 

VioScroUUp(0, 0.-1, &(char)Cell,0); 
puts(”Press Both Mouse Buttons To Exit"); 

MouD rauPt r ( MouHand I e ) ; 
do C 

sprintf(OutStr,"X*X2d Y*X2d", MouEvent. col, MouEvent. row) ; 
VioW^tCha^St^<CXJtSt^,strlen(OutStr),0,0,0); 

MouR eadE ven t Que ( &HouE vent , &Wa i t Opt i on , MouH and I e ) ; 

> whi I e( (MouEvent .fs A 0x14) 1= 0x14); 

MouC I ose ( MouHand I e) ; 
puts ("Have a Mice Dayf"); 
au_end(); 

> 
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/* File SAVE.C V 

#include <stdio.h> 

#include <8tdlib.h> 

#include <strfng.h> 

#include ••util.h” 

typedef int UORD; 

typedef long DUORO; 

void mainCint argc,char *argvCJ); 

void main(argc,argv) 
int argc; 
char *argv[]; 

< 

char f name [130]; 
int ha; 

if (argcss2) strcpy(fname,argv[1] ); 
else { 

printf (“Enter file to save: "); 
scanf ( “Xs” , f name) ; 

> 

au_init(&ha); 

printf (“Saving Xs. . .\n“,fname); 
au_save( ha , f name) ; 
au_out put ( ha , 0x8000 ) ; 
au_end(ha); 

> 
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/* File ZOOH.C V 

#include <stdio.h> 

#include <conio.h> 

#include "util.h** 

unsigned char huge ipict512] [512] ; 
unsigned char huge op1c[512] [512]; 
void maln(vold); 
void zoom(lnt); 

Int ha; 

vold maInO 
C 

Int 1; 

/* Initialization V 
au_1n1t(&ha); 
for (1=0;1<3;1 -m-) i 

prIntfC "Zooming buffer Xd...\n",1); 
zoom(l); 

> 

au_end(ha); 

> 

void zoomCbuff) 

Int buff; 

< 

Int 1J,x,y; 

au_readscr(ha, 1p1c,buff ); 
for <1=0;1<512;i^»2) 

for (j=0; j<512;j>=2) { 
x*128+1/2; 
y=128^j/2; 

opIcCIlCJ] * ip1c[x][y]; 

opIcCll [j>11* dpic[x] [y]+ip1c[x] [y+1] )/2; 

opic[1+1] [jl= (1p1c[x] [y]+1p1c[x+13 [y])/2; 

opic[1*»-1] [j+1]= dpiclx] ty] + 1pic[x+1l [y^1])/2; 

> 

au_wr1te8cr(ha,op1c,buff); 

> 
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